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The  mechanisms,  kinetics,  and  reversibility  of  ortho- 
phosphate adsorption  on  synthetic  goethite  and  soils  were 
investigated.  Goethite  was  prepared  by  ageing  of  the  precip- 
iPsTs  which  appears  after  mixing  PeCl2  and  NaOH  solutions. 
Both  the  increase  in  Fe/OH  ratio  in  the  suspension  and  in 
ageing  temperature  favored  a higher  degree  of  goethite  crys- 
tallization within  a short  period  of  time.  With  either 
OH/Fe  = 6 in  the  suspension  at  room  temperature  or  OH/Fe  = 3 
at  55°C  and  ageing  for  one  week,  good  goethite  yield  was 
obtained.  Infrared  bands  characteristic  of  goethite  were  at 
3200  cm  ^ (OH  streching)  and  at  890  cm~^  and  790  cm~^  (both 
are  Fe-OH  bending  vibration).  The  goethite  structure  was 
confirmed  by  X-ray  diffraction  intensities  of  the  4.19  8 and 
2.70  8 peaks.  The  presence  of  phosphate  at  the  beginning  of 
goethite  ageing  weakened  the  goethite  structure.  I-Then  F/Fe  = 
3.2  (regardless  of  OH/Fe  ratio)  the  bond  (Fe)-  0-P  vibration 
at  1000  cm  ^ was  predominant,  thus  preventing  the  formation 
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of  Fe-OPI  bond  even  after  ageing  of  a suspension.  When  the 
addition  of  phosphate  was  made  at  the  end  of  the  goethite 
ageing,  absorption  of  phosphate  was  determined  as  surface 
binding  through  binuclear  bridging  of  the  HPO^^  ion  which 
was  identified  by  the  presence  of  bands  at  1120,  1085,  and 
1030  cm  The  phosphate  was  assumed  to  penetrate  the 
goethite  structure  whenever  the  vibrational  band  at  1000  cm~^ 
was  present. 

Using  synthetic  goethite  as  supporting  medium,  it  was 
found  that  the  Langmuir  and  Freundlich  equations  could  be 
used  to  describe  the  relation  between  the  amount  of  P adsorbed 
and  that  remaining  in  the  equilibrium  solution  for  particular 
periods  of  reaction  time. 

A kinetic  model  was  proposed  to  describe  the  change 
with  time  of  each  of  the  following  ion  forms:  1)  the  free 
ions  in  solution,  2)  the  physically  adsorbed,  3)  the  revers- 
ible chemically  adsorbed  ions,  and  4)  the  irreversible 
chemically  adsorbed  ions . The  amount  of  phosphate  adsorbed 
increased  with  the  increase  in  initial  P concentration  and 
the  adsorption  reaction  time.  The  time  required  for  the 
goethite  phosphate  systemi  to  reach  an  equilibrium  state 
increased  with  the  initial  P concentration.  Generally,  the 
equilibrium  state  was  reached  after  18  hours  for  goethite- 
solution  and  22  hours  for  soil-solution  systems.  The  adsorp- 
tion on  goethite  was  increased  by  both  multivalent  cations 
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and  the  concentration  of  the  supporting  electrolyte.  The 
potential  binding  energy  constant  increased  (from  0.70  to 
4.60  ml/  yg  P)  as  the  reaction  time  increased,  but  decreased 
as  the  cation  valence  of  the  electrolyte  was  increased. 

At  the  same  initial  P concentration,  the  amount  of  P 
adsorbed  decreased  almost  linearly  with  change  in  pH  accord- 
ing to  the  relation  yg  P/g  goethite  = -446  pH  + 5680.  The 
amount  of  P desorbed  over  a wide  range  pPI  or  supporting 
electrolyte  remained  nearly  constant  when  the  desorbing  time 
was  greater  than  or  equal  to  6 hours.  The  P desorbed  from 
goethite  and/or  soils  increased  both  with  the  initial  equi- 
librium P concentration  employed  and  the  reaction  time.  At 
the  same  initial  P adsorbed  on  goethite,  the  amount  of  P 
desorbed  from  goethite  decreased  when  the  pH  was  below  5.5 
and  increased  when  the  pH  was  greater  than  6 . Phosphate 
adsorption  on  goethite  was  found  to  induce  a net  increase  of 
negative  charge  so  that  the  zero  point  of  charge  declined 
from  pH  5.8  to  5.2. 
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INTRODUCTION 


Increase  in  urbanization  makes  it  necessary  to  increase 
both  the  quantity  and  the  quality  of  agricultural  production 
in  the  tropical  region.  The  implementation  of  an  adequate 
agricultural  policy  involves  improvements  of  the  actual 
level  of  technology.  This  implies  a clear  understanding 
that  the  macroscopic  relation  be  sought  between  the  major 
factors  producing  food  for  human  population  through  proper 
managements  of  soils,  plants,  animals,  and  insects  under 
diverse  climates  and  social  activities.  Modeling  of  the 
agricultural  system  is  a useful  tool  for  understanding  a 
united  approach  for  all  major  factors  governing  the  system. 
With  a suitable  model,  it  should  be  possible  to  envisage 
needed  change  in  a particular  factor  in  order  to  give  results 
as  close  as  possible  to  a reasonable  expectation.  Only  some 
microscopic  relationships  within  the  soils  will  retain  our 
attention  in  this  study.  In  a soil-solution  system,  the 
dominant  soil  phenomena  taking  place  simultaneously  are  mass 
transport,  diffusior^  adsorption-desorption , precipitation- 
dissolution,  and  microbial  immobilization  and  mineralization. 
The  dominance  of  each  soil  phenomenon  depends  on  the  soil 
structure  and  texture,  organic  m.atter  content,  water  content, 
temperature,  and  the  ionic  suite  (types  and  concentrations). 
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Upon  fertilization,  tropical  soils  may  not  react  in  the 
same  way  as  do  the  temperate  soils  because  of  their  different 
mineralogical , and  chemical  properties.  Among  the  major 
elements,  phosphorus  is,  after  nitrogen,  the  most  deficient 
nutrient  and  its  availability  is  strongly  dependent  on  the 
mineralogical  composition  of  the  soil  (type  of  clay,  and 
metal  oxides).  The  low  availability  of  phosphate  to  tropical 
plants  is  due  both  to  the  high  soil  phosphate  fixing  capacity 
and  to  precipitation  of  fertilizer  phosphorus  from  the  soil 
solution.  Not  all  forms  of  phosphorus  bound  on  the  oxide 
surfaces  are  held  with  the  same  degree  of  strength.  V/ith 
time,  as  phosphorus  is  depleted  from  the  solution,  some 
phosphate  may  be  replenished  by  release  from  the  solid 
phase.  The  magnitude  of  this  release,  vjith  respect  to  time, 
is  dependent  on  soil  characteristics.  The  initial  concen- 
tration in  solution  is  not  a sufficient  measure  of  phosphorus 
availability.  In  tropical  soils,  iron  oxides  are  quite 
important  in  determinig  the  solid  phase  capacity  to  fix  and 
supply  phosphorus  to  the  soil  solution.  How  some  factors, 
such  as  OH/Fe  ratio,  temperature,  and  phosphate  concentra- 
tion affect  the  formation  of  goethite  are  herein  investi- 
gated. Infrared  spectroscopy  will  he  used  to  identify  the 
effects  of  soluble  phosphate  on  goethite  structure  and  the 
nature  of  phosphate  binding  on  goethite.  Synthetic  goethite 
as  well  as  three  soils  will  be  used  to  investigate  some 
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factors  affecting  the  time  dependence  and  the  degree  of 
reversibility  of  orthophosphate  adsorption. 


CHAPTER  I 

REVIEV/  OE  SOME  MODELS  DESCRIBING  EACTORS 
AFFECTING  NUTRIENT  AVAILABILITY 

Nutrient  Potential 

The  concept  that  availability  of  a nutrient  can  be 
described  by  its  potential  activity  was  first  brought  to 
general  notice  by  Schofield  (1955).  His  concept  was  that  in 
comparing  two  soils  of  different  water  holding  capacity 
which  contain  the  same  amount  of  available  water,  the  one 
having  the  lower  capacity  has  a lower  potential  and  the 
water  is  easier  to  extract  than  from  the  soil  of  higher 
holding  capacity.  Their  pF  value  defines  the  energy  v;ith 
which  water  is  held  on  the  soil  particle  surface.  The 
potential  given  by  the  pE  value  provides  the  basis  for 
quantitative  evaluation  of  the  availability  of  water.  As 
water  is  taken  out  of  a soil,  it  is  replaced  by  lateral 
movement  of  groundwater.  Similarly,  as  a nutrient  is  taken 
from  the  soil  solution,  it  is  replaced  by  other  ions,  by 
desorption,  or  by  diffusion  as  well  as  mass  flow.  Eventually, 
it  is  not  the  amount  of  a nutrient  in  a soil  that  primarily 
controls  the  uptake  of  that  nutrient  by  the  plant,  but  the 
work  required  to  withdraw  it  from  the  solution.  This  work 
may  be  related  to  Gibb's  free  energy  (G^)  if  the  uptake  is 
mainly  from  the  soil  solution.  To  derive  an  expression  for 
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the  nutrient  potential,  three  postulates  have  been  advanced 
by  different  workers. 

Postulate  1 

Where  the  solution  is  at  constant  temperature  and 
pressure,  and  when  the  system  is  at  equilibrium,  is 
uniform  throughout  the  soil  solution  so  that 

G.  = G?  + RTLn  a^  + zF'l'f  (1) 

where  G?  is  the  standard  molar  free  energy  of  an  ion  i in 
the  solution  relative  to  an  arbitrarily  zero  of  the  electric 
potential,  a^  is  the  activity  of  this  ion  in  solution,  z is 
the  valence,  F is  the  Faraday  constant,  and  is  the 
electric  field  effects. 

Postulate  2 

When  a soil  is  in  equilibrium  with  a solution,  the 
electro— chemical  potential  of  the  ion  is  constant  throughout 
the  system  (soil-solution).  Then,  the  partial  molar  free 
energy  of  any  ion  in  the  soil  complex  can  be  determined  by 
analyzing  the  solution.  This  free  energy  of  the  ion  in 
solution  is  assumed  to  be  constant  even  if  the  solution  is 
separated  from  the  solid  phase  by  centrifugation.  As  a 
result,  if  the  ion  is  taken  out  of  the  field  influence,  then 
zF4'  = 0 where  z,  F,  and  4'  are  as  described  in  the  previous 
section.  Then  the  relationship  in  Eq . (1)  can  be  reduced  to 


6 


G.  = G?  + RTLn  af  (2) 

Postulate  3 

This  postulate  is  that  potentials  are  determined  by 
transfers  of  ions.  In  order  to  maintain  the  electroneu- 
trality, equal  transfer  must  occur  for  quantities  of  posi- 
tive and  negative  ions  in  one  direction  or,  alternatively, 
transfer  of  one  ion  in  the  first  direction  occurs  as  an 
equivalent  quantity  of  ions  of  the  same  sign  moves  in  the 
opposite  direction  (Barrow  et  al.  1965).  It  is  also 
assumed  that  divalent  cations  (Ca  and  Mg)  dominate  m.ost  soil 
surfaces  as  exchangeable  cations  and  in  the  solution,  except 
for  saline  soils.  For  a situation  in  which  transfer  of  Ca^"*" 
is  in  one  direction  and  K’*’  in  the  other,  the  net  change  in 
free  energy  is 

^^Ca,K  " ^^Ca,K  ^Ca  ^ 

For  phosphorus,  over  the  range  of  pH  v;hich  exists  in 
soil,  two  forms  of  the  orthophosphate  ions  are  in  equilib- 
rium; thus 

H PO"  + H O ^ = = = = = = i HPO^^~  + H.O'^  (4) 

pK=7.2  ^ 

Assuming  that  plants  absorb  m.ainly  H^PO^  ion  and  not  HPO^~, 
then  a cation  such  as  Ca  must  accompany  the  H^PO^.  The 
equation  for  the  transfer  of  both  ions  from  the  equilibrium 
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^^Ca(H2PO^)2  ^ RT(l/2pCa  + pH^PO^^)  (5) 

A common  use  of  the  relationship  is  the  Schofield's  phospho- 
rus potential,  where  I is  given  by  I = 1/2  Ca^"*"  + pH^PO^. 

Capacity  and  Intensity  Relationship 
Phosphorus  ions  pass  from  one  phase  to  another  as  a 
result  of  chemical  potential  differences  between  the  phases. 
In  general,  any  substance  tends  to  pass  from  a region  of 
higher  chemical  potential  to  a lower  one.  The  ability  of 
the  solid  phase  to  supply  phosphorus  to  the  solution  phase, 
as  it  is  depleted,  can  be  termed  the  capacity.  For  any  ion 
such  as  phosphorus,  at  equilibrium,  the  phosphorus  buffering 
capacity  (PBC)  is  a soil  characteristic  where  PBC  is  dAP/dl. 
The  term  with  I = l/2pCa^'^  + pH^PO"  and  AP  is  the  gain  or 
loss  in  phosphorus  by  the  soil.  The  relationship  between 
capacity  and  intensity  factors  is  expressed  as  Q/I  curves 
which  are  generally  composed  of  two  parts  with  linear  and 
curvilinear  portions.  In  the  Q/I  relationship  for  potassium, 
the  curved  portion  is  attributed  to  a number  of  sites  which 
have  specific  affinity  for  potassium  at  lov;  concentration 
and  the  linear  part  is  associated  to  non-specific  sites, 
Beckett  (1971).  Ke  also  believed  that,  in  general,  the 
curved  part  of  the  Q/I  relation  can  be  represented  by  a 
Langmuir  adsorption  isotherm  whereas  the  linear  part  fits 
the  Gapon  relationship. 
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Becketrlr  (1971)  affirTned  "that  "the  non-linear  part  of  the 
curve  represents  a certain  region  where  there  is  a definite 
limited  number  of  sites  on  the  exchangeable  surfaces  and 
exhibits  a selective  binding  power  for  the  adsorbing  ion: 

The  Langmuir  adsorption  isotherm  could  be  used  to  describe 
such  a situation.  In  addition,  if  Eq.  (3)  is  recalled,  the 
activity  ratio  (AR)  is  K'^/CCa)^^^  of  the  ions  in  solution 
which  is  the  product  of  the  concentration  of  the  exchange- 
able ions  and  the  Gapon  constant  (k^). 

Energy  of  Adsorption 

Partition  function 

To  be  effective,  the  collision  between  molecules  and 
the  collision  surface  must  provide  a certain  m.inimum  amount 
of  energy  called  the  activation  energy.  Since  the  activated 
complex  is  a transitory  species,  the  equilibrium  constant 
cannot  be  measured  experimentally.  However,  the  partition 
function  arises  from  the  quantum  theory  that  a molecule  can 
exist  only  in  states  with  definite  energy  limits.  from 
Boltzman  distribution,  it  is  recognized  that 

where  is  the  total  number  of  molecules,  N . is  the  number 

o/\ 

of  molecules  at  zero  energy,  E^  is  the  activation  energy  at 
the  state  i , and  g^  is  a constant.  In  this  case,  also 
quantity  Q = Zg^^^exp ( -E^^/kT ) which  is  the  partition 
function. 
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Electric  Field  Effect 

The  following  treatment  was  first  derived  for  solid-gas 
interaction  and  is  known  as  Thin's  theory  reported  by 
Brunauer  (1943).  The  strength  of  the  electric  field  sur- 
rounding the  adsorbent  is  given  by  E — 9 where  e^^  is 

the  charge  distribution  of  the  absorbent,  e is  the  dielec- 
tric constant  of  the  gas,  and  r is  the  distance  from  the 
surface . 

The  force  (F)  acting  on  an  induced  dipole  is 
F = ke^v/4  nr^^"^^ 

where  v is  a constant  and  k = e-1/ e . The  adsorption  poten- 
Tis-l  is  the  force  multiplied  by  the  distance  through  which 
its  acts  (F.Ar). 

Heat  of  Adsorption 

Assuming  no  work  due  to  phase  change  is  done  during  the 
adsorption-desorption,  the  molar  energy  of  the  ion  in  free 
solution  (u^)  and  the  energy  of  the  ion  in  the  adsorbed 
phase  (Ug)  are  related  to  the  loss  of  ion  from  the  solution 
due  to  adsorption  as  shov/n  by 
AH  = b(u^  - 

vjhere  b is  the  number  of  miole  adsorbed,  and  AH  is  the  integ- 
ral heat  of  adsorption. 

Differential  Heat  of  Adsorption 

Assuming  that  the  adsorption  process  is  reversible  AH 
can  be  obtained  from  data  at  two  temperatures  (T^  and  T^ ) by 
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using  the  Clausius-Clapeyron  equation  for  a constant  surface 
coverage,  namely 


^2 

AH 

1 1 

R 

1 

1 

ro  1 

where  a^  and  a^  are  equilibrium  concentrations  at  tempera- 
ture and  respectively,  and  R is  the  molar  gas  constant 

Some  Adsorption  Isotherm  Models 
Freundlich  Equation 

The  Freundlich  equation  was  first  introduced  in  an 
empirical  form  (Bach  and  V/illiams  1971).  It  assumed  that 
the  energy  of  adsorption  decreases  exponentially  with  increas 
ing  saturation  of  the  surface.  The  equation  is 

b/m  = ka^^^  (9) 

where  b/m  is  the  amount  of  phosphorus  adsorbed  per  unit 
weight  of  soil  at  the  equilibrium  concentration  termed  a, 
and  where  k and  n are  constants.  The  logarithmic  form  of 
Eq.  (9)  is 

log(b/m)  = log  k + 1/n  log  a (10) 

The  plot  of  log(b/m)  versus  log  a should  yield  a 
straight  line.  This  equation  is  valid  only  in  a limited 
range  of  concentrations.  By  taking  into  account  the  ini- 
tially exchangeable  phosphorus  and  the  ability  of  the  soil 
to  lose  or  gain  phosphorus  (AP)  during  the  equilibrium 
reaction,  the  following  modification  was  introduced 

AP  = Aa  = ka^'^’^-  e (11) 
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where  Aa  is  the  amount  of  phosphorus  gained  or  lost,  e is 
the  phosphorus  initially  present,  and  n is  a constant. 

Modified  Freundlich  Equation 

In  a colloid-solution  system,  the  adsorption  is  time- 

dependent.  Kuo  and  Lotse  (1972)  introduced  a time  factor 
* }c  t 

1-e  2 after  assuming  that  the  constant  is  small.  They 

found  that 

b/m  = ka^t^'^^  (12) 

where  a^  is  the  initial  phosphorus  concentration.  Their 
study  indicated  that  the  rate  constant  increased  with  an 
increase  in  concentration. 

Langmuir  Equation 

The  derivation  of  the  Langmuir  equation  is  based  on 
the  following  assumptions:  1)  the  energy  of  adsorption  is 

constant  and  independent  of  the  degree  of  coverage,  2)  there 
is  not  interaction  between  adjacent  adsorbed  molecules  on  an 
homogenous  surface.  If  the  system  is  in  dynamic  equilib- 
rium, this  results  when  the  rate  of  adsorption  is  equal  to 
the  rate  of  molecules  escaping  from  sorption  surface,  so  that 

^1®  ■ ^2  ^13) 

where  and  k^  are  the  rate  constants  of  adsorption  and 
desorption,  P is  the  gas  pressure,  and  6 is  the  fraction  of 
the  surface  coverage  by  the  gas.  By  rearrangement,  the 
relationship  is 
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(14a) 


and  by  analogy 


for  soil-solution  system,  where  a is  the  ion  concentration 
in  solution.  The  above  equation  can  be  put  in  the  linear 
form  as  follows 


where  k is  the  adsorption  energy  constant  and  s is  the 
adsorption  maximum.  Different  workers  have  observed  that 
there  is  not  a linearity  between  a/b/m  and  a.  This  may  be 
caused  by  at  least  two  types  of  sites  of  adsorption  having 
different  energies  of  adsorption.  To  fit  this  criterion, 
the  Langmuir  equation  was  modified  to  be 


where  k^,  s^  are  constants  for  region  1,  and  k^,  S2  are 
constants  for  region  2 (Syers  et  al.  1973). 

Effect  of  Surface  Heterogeneity  on  Adsorption-Desorption 
Langmuir  and  Freundlich  equations  can  be  obtained  from 
Toth's  equation  (Jossen  et  al.  1978),  for  homogeneous 
surface.  Toth's  equation  is  obtained  by  integrating  the 
equation 


a/b/m  = a/s  + 1/ks 


(14c) 


(14d) 


(15) 


Where  6=0,  the  Freundlich  equation  is  obtained.  Where 
6=1,  the  Langmuir  type  of  equation  is 
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(16) 


where  6 = d exp(-E/RT),  a is  the  equilibrium  concentration, 


Another  factor  is  that  the  energy  of  activation  (E)  is 
constant  for  a homogeneous  surface  but  varies  on  a hetero- 
geneous surface  with  the  degree  of  coverage  (Jossen  et  al. 
1978).  Assuming  that  the  free  energy  of  activation  for 
adsorption  varies  linearly  with  coverage  of  the  surface, 
Lingstom  et  al.  (1970)  proposed  the  following  model 


where  A is  the  adsorbate,  S is  the  surface,  and  and  k2 
are  the  rate  constants.  They  deduced  that  the  rate  of 
adsorption  process  is 


de/dt  = k^(l-0)(l-l/26)e  + k2[(l-l/20)e  ®^-l/20^e^®] 

(18) 


Another  useful  relationship  is  the  Elovich  equation 
where  the  activation  energy  is  a linear  function  of  the 
amount  adsorbed 


Aharoni  and  Ungarish  (1977)  modified  this  equation  by  intro- 
ducing the  fact  that  the  heterogeneous  surface  is  comprised 
of  a large  number  of  homogeneous  regions  having  unequal 
number  of  adsorption  sites. 


b the  amount  adsorbed  and  b is  the  maximum  adsorption. 


(17) 


where  0 = b/b°°. 


E = E + ab 
o 


(19a) 
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= Eq  + RTln(gy^/n^  + y)  (19b) 

where  = db/dE,  g and  y ai'e  constant,  and  n^  is  the  number 

of  sites  of  energy  , where  E^  is  the  activation  energy 

characteristic  of  a region.  They  also  assumed  that,  at  any 

moment,  adsorption  takes  place  preferentially  on  the  region 

that  has  the  lowest  activation  energy  at  that  moment. 

Considering  the  case  where  equilibrium  is  attained  at  a 

region  when  y^  = y^^j  then  it  follows  that  y^^  = k^exp(E/RT). 

The  overall  rate  of  adsorption  is  given  by 

db/dt  = k N exp(-E./RT)  and  N.  = /n„dE  (20) 

o t ^ t t E 

where  is  the  number  of  sites  with  an  energy  E^  at  time  t 
and  n^  is  the  number  of  sites  with  an  activation  energy 
between  E and  E + dE. 

Generalization  of  Adsorption  Isotherms 
Freundlich  and  Langmuir  equations  can  be  extended  for 
cases  of  competitive  adsorption  (Jaroniec  and  Toth  1976; 
Digiand  et  al.  1978).  In  this  case,  the  Freundlich  type  is 
given  by 

b^  = k(Za^)  (21) 

i 

and  the  Langmuir  type  by 

b^  = ka^/ (k^  + k^a^  + ^2^2  (22) 

which  is  a partial  isotherm  of  a^  relative  to  the  total  ions 
present.  A binary  system  can  be  simplified  if  the  following 
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assumptions  are  made,  1)  the  surface  is  formed  by  a collec- 
regions  each  being  characterized  by  3.  binding  energy 
(E),  2)  the  number  of  surface  atoms  Ih  with  energy  E follow 
the  Boltzman  distribution,  which  is  described  by 

Ni  = N^exp(-e/RT)  (23a) 

and  3)  all  adsorbed  molecules  stay  on  the  surface  for  the 
same  average  length  of  time  6 so  that 

6 = 6^exp(E/RT)  (23b) 

During  the  time  6,  the  fluctuation  of  the  number  of  adsorbed 
and  desorbed  molecules  compensate  between  themselves  randomly 
(Vlad  and  Segal  1979).  They  considered  that  the  adsorption 
energy  is  an  increasing  function  of  the  extra  energy  e. 

From  Mclaurin  development,  the  extra  energy  can  be  expressed 
as 

e = Za^(E  - E )^  (24) 

where  E^  is  the  smallest  adsorption  energy  corresponding  to 
the  value  of  zero  energy  and  n is  generally  integers  of  1 or 

2. 

The  energy  distribution  X(e)  is  expressed  by  Van  Dongen 
approximation  as 

X(e)  = exp(Z  a^e)  i = 0,  1,  or  2.  (25a) 

i 

or  generally  as 

N , N 

X(e)  = 1/kT  I na  (E  - E )^~-^exp[-l/kTZa  (E  - E )^] 
n = T ^ ^ m 


(25b) 
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The  general  isotherm  is  expressed  as 

b(a,T)  = /^b^(a,e)  X(s)  de  (26) 

where  b(a,T)  is  the  fraction  of  the  total  surface  covered  at 
(a,e),  b^(a,e)  is  the  local  isotherm  which  may  be  analogous 
to  an  isotherm  equation  for  an  homogeneous  surface  or  region, 
is  the  range  of  possible  variation  of  adsorption  energy 
assuming  that  = [0,«>]. 

Temperature  Dependence  of  the  Rate  Constant 

According  to  the  Arrhenius  equation,  it  is  known  that 

^ads  " Aexp(-E^/RT)  (27a) 

where  is  the  energy  of  activation,  R is  the  gas  constant, 
and  T is  the  absolute  temperature. 

Erom  the  transition  state  theory,  this  relationship  is 
extended  to 

kads  = (kT/h)exp(AS^/R)exp(-AH^/RT)  (27b) 

= (kTZh)exp(AG  /RT)  (27c) 

where  AS  is  the  entropy  of  activation,  AH  = E - RT  which  is 
the  enthalpy  of  activation,  and  is  the  Gibb's  free  energy 

of  activation. 

Soil  Phosphorus  Reaction  Mechanisms 
Plant  responses  to  soil  phosphorus  are  a function  of 
the  solubility  of  phosphorus.  Any  factor  altering  the 
solubility  of  phosphorus  will  also  alter  the  plant  response 


17 


(Hemwall  1957).  The  solubility  of  phosphorus  depends  on 
several  factors  such  as  the  composition  of  the  soil  solu- 
tion, the  pH  of  the  solution,  and  temperature  of  the  soil- 
solution  system. 

Calcareous  and  Neutral  Soils 

The  adsorption  of  phosphorus  on  calcareous  surfaces  can 
take  place  by  replacement  of  water  molecules,  bicarbonate, 
and  certain  other  anions  or  cations.  The  relative  strength 
of  the  phosphate  ion  adsorption  depends  on  the  solubility  of 
the  compound  at  the  calcium  surface  (Kuo  and  Lotse  1972). 

The  precipitation  of  phosphorus  may  be  due  to  the 
formation  of  a whole  series  of  insoluble  calcium  phosphates. 
Some  of  these  with  associated  solubility  products  expressed 
as  the  pK  value  are  reported  by  Lindsay  and  Moreno  (1960). 


Compounds 

Chemical  formula 

pK 

Calcium  phosphate  anhydride 

CaHPO^^ 

6.66 

Calcium  phosphate  dihydrate 

CaHPO^^.  2H2O 

6 . 56 

Octocalcium  phosphate 

Ca^^H(P0^^)3.  SH^O 

46.91 

Hydroxyapatite 

Caio(PO^^)6.(OH)2 

113 .70 

Fluorapatite 

Caio(PO^^)6.p2 

118.40 

In  calcareous  soil,  hydroxyapatite  and  fluorapatite  are 
the  major  phosphorus  compounds,  whereas  in  neutral  soil, 
octocalcium  phosphate  becomes  important  and  in  moderately 
acidic  soils  dicalcium  phosphate  may  occur. 
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Acid  Soils 

Iron,  aluminum  and  pH  are  the  main  factors  controlling 
phosphorus  solubilities;  In  acid  soils,  the  solubility  of 
phosphorus  increases  with  decreasing  free  iron  and  aluminum 
activities  and  with  increasing  pH.  The  form  of  phosphate 
sorbed  by  acidic  soils  was  recovered  (90%)  as  iron  and 
aluminum  phosphates  (Ghani  and  Islam,  1946).  Yuan  et  al. 
(1960)  reported  that  up  to  80%  of  the  added  phosphorus  to 
acid  sandy  soils  was  present  as  aluminum  phosphate  and  10% 
as  iron  phosphate,  but  when  the  reaction  temperature  was 
increased  more  iron  phosphate  was  formed. 

After  reaction  of  phosphorus  with  soluble  Al^'*',  a 
microscopic  examination  showed  an  hexagonally  shaped  crystal 
in  which  the  interplanar  spacing  was  similar  to  those  of 
palmerite  (Haseman  et  al.  1950).  In  general,  some  of  the 
main  forms  in  which  phosphorus  can  precipitate  with  iron  and 
aluminum  are  given  below. 


Compounds 

Chemical  formula 

pK 

30.5 

Variscite 

AlPO^^.  2H2O 

Strengite 

FePO^^.  2H2O 

34.3 

Taranakite 

(K,NH^^)3HgAl^(PO^^)g.l8H2  0 

176 . 0 

Palmerite 

HK2Al2(PO(^)g.  7H^0 

Since  Al  and  Fe  atoms  are  part  of  the  surface  colloids, 
they  react  with  phosphates.  Whether  the  process  is 
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precipitation  or  adsorption  depends  on  the  size  of  the 
metal  polymer  and  the  pH  of  the  phosphate,  and  its  concentra 
tion.  In  a moderately  acidic  medium,  with  a high  phosphorus 
concentration,  the  reaction  process  may  be  typically  precipi 
tation,  resembling  that  reported  by  Hsu  (1965), 

[Alg(0H)^2^^^  + 6 H^PO”  Alg(OH)^2^H2PO^^)g  (28) 

Specific  Adsorption 

General  Description 

In  a solid-solution  system,  adsorption  of  molecules  (or 
ions)  occurs  when  there  is  a change  in  phase  from  the  free 
state  (in  solution)  to  the  bound  state  at  the  interface. 

Ions  deposited  at  the  surface  likely  orient  to  form  the 
Stern  layer  as  some  ions  may  approach  closely  to  the  surface 
structure.  In  this  case,  the  ions  are  said  to  be  specifi- 
cally adsorbed  on  the  surface. 

Non-specifically  adsorbed  ions  are  either  in  the  dif- 
fuse Gouy-Chapman  region  separated  from  the  solid  surface  by 
at  least  one  molecule,  or  electrostatically  bound  to  that 
surface.  Specifically  adsorbed  ions  are  in  the  coordination 
shell  of  the  surface  atoms  and  are  maintained  there  through 
chemical  binding  (covalent  or  coordinate  binding).  Since 
specific  adsorption  of  cations  or  anions  occurs  even  when 
the  surface  possesses  a net  positive  or  negative  charge, 
respectively,  there  must  be  an  electrostatic  contribution 
due  to  polarization  of  the  ion  or  molecule. 
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Both  chemical  and  electrostatic  attractions  contribute 
to  the  energy  of  adsorption.  The  magnitute  of  this  energy 
of  adsorption  determines  the  degree  of  reversibility  and  the 
amount  adsorbed.  Since  the  reaction  is  pPI  dependent,  at  any 
pH,  there  is  a maximum  adsorption  of  anions  and  when  these 
maxima  are  plotted  against  pH,  the  curve  is  termed  the 
adsorption  envelope  which  is  described  by  Kingston  et  al. 
(1967)  as 


where  A is  the  amount  of  the  ion  adsorbed  per  unit  weight  of 
adsorption,  a is  the  degree  of  dissocation  of  acid  anion,  K 
is  the  dissociation  constant  of  the  most  highly  charged 
anion  that  is  adsorbed,  and  V is  the  amount  of  ion  adsorbed 
at  the  maximum  level. 

Why  som.e  of  the  adsorbed  anions  such  as  phosphate  may 
be  irreversibly  held  is  explained  by  certain  postulated 
mechanisms  of  phosphate  adsorption.  The  removal  of  phos- 
phate from  free  solution  is  assumed  to  proceed  through  the 
replacement  of  coordinated  H^O  groups  and/or  some  OH”  ions. 


A = 2Vn/U  - a)  = 2V 


Kun 


(29) 


(K  + [H])^ 


OH 


OH 


0-P  = 0 


M 


OH  + OH 


(30) 


OH 


OH 


OH 


where  M is  a metal  (Fe  or  Al ) , and  bonded  by  a coordination 
link . 


21 


Infra-red  Study  of  Phosphate  Specific  Adsorption 

Theory  of  Infra-Red  Spectroscopy 

If  a molecule  is  placed  in  an  electromagnetic  field,  a 
transfer  of  energy  from  the  field  to  the  molecule  occurs 
when  Bohr's  frequency  condition  is  satisfied:  AE  = hv, 

where  AE  is  the  difference  in  energy  between  two  quantized 
states,  h is  Plank's  constant,  and  v is  the  frequency  of  the 
light  wave.  Pure  rotational  vibrations  are  usually  observed 
in  the  microwave  vibrational  spectra  in  the  infrared  whereas 
electronic  spectra  are  in  the  visible  and  ultraviolet.  The 
infrared  spectra  originate  in  the  transition  state  existing 
between  two  vibrational  levels  of  the  molecules.  From  a 
classical  point  of  view,  a vibration  is  active  in  the  infra- 
red spectrum  if  the  dipole  moment  of  the  molecule  is  changed 
during  the  vibration.  The  dipole  moment  Q is  related  to  the 
strength  of  the  electric  field  by  Q = aE  and  Q is  a vector 
whose  direction  is  the  line  between  the  center  of  gravity  of 
the  protons  and  electrons . Let  a diatomic  molecule  be 
represented  by  two  masses  m^  and  m^  moving  along  the  mole- 
cule's axis  with  displacements  of  and  respectively. 

T'he  displacements  of  the  two  atoms  are  induced  by  forces 
which  can  be  obtained  through  Hooke's  and  Newton's  law. 


KCX^  - X^)  = m^ 


dt 


dt 


2 

2 


(31) 
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The  solutions  of  these  equations  of  motion  are: 

= Acos(27TVt  + a)  and  - Acos(2iTVt  + a)  (31b) 

After  differentiating  and  substituting  back  in  Eq . (31),  it 
can  be  found  when  solving  for  v (Nakamoto  1978;  Colthup 
et  al.  1975) 


(32) 


where  k is  the  force  constant  and  v is  the  frequency  of 
vibration. 

The  experimental  observations  of  frequency  of  vibration 
of  a crystal  are  due  to  bond  vibrations  of  atoms  within  the 
unit  cell.  Such  vibration  m.ay  help  us  to  determine  the 
nature  of  the  atoms  composing  the  unit  cell  by  comparison  to 
known  polyatomic  vibrations.  Because  of  the  interactions, 
the  symmetry  of  a molecule  is  generally  lower  in  the  crystal- 

. q 

line  state  than  in  the  isolated  state.  The  isolated  PO^ 
ion  is  of  tetrahedral  form  (Nakamoto  1978): 


0 


/ 


t 

0 


The  above  tetrahedral  structures  predict  two  infrared  active 
fundamentals,  one  is  stretching  (v^)  and  the  other  (^^)  is 
bending . 
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Infrared  Identification  of  the  Forms  of  P Adsorption 

From  their  studies,  Kingston  et  al.  (1974)  proposed  the 
reaction 


'^0H2 


-2 


0 


(33) 


Reversibly  adsorbed  P 


Irreversibly  adsorbed  P 


Parfitt  et  al.  (1975)  identified  a form  of  phosphorus 
which  can  be  considered  to  be  completely  within  the  goethite 
structure,  since  the  phosphate  is  bound  with  two  Fe  atoms 
and  is  H-bonded  to  a third  Fe. 


Fe 


0 HO 


0 


0 - Fe 


0 - Fe 


Fe 


OH 0 


or 


> 


0 


.0  - Fe 


'0  - Fe 


(34) 


Such  a chemical  binding  of  phosphate  on  the  goethite  surface 
is  said  to  be  specifically  adsorbed,  and  can  be  considered 
as  the  formation  of  a new  solid  phase  or  growth  of  the  solid 
phase . 

Colloid  or  Soil  Surface  Effects 

Substances  dissolved  in  soil  solution  can  move  by 
molecular  or  ionic  diffusion  resulting  from  a concentration 
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gradient  within  the  solution,  or  by  mass  flow  of  the  soil 
solution.  This  is  complemented  by  sorption  of  the  ions  onto 
the  soil  surface,  or  by  a combination  of  these  factors.  The 
above  phenomena  occur  simultaneously  and  are  described  by 
the  following  equation  of  de  Camargo  et  al.  (1979): 


dt  ^^2  dx  dt  dx 


(35) 


where  a is  the  ion  concentration  in  the  soil  solution,  b is 
the  amount  of  the  ion  adsorbed,  x is  the  distance  of  the  ion 
from  the  electric  potential  effect  due  to  the  surface  charge 
on  the  charged  solute.  In  Eq . (35),  the  constants  are  A 
which  is  a coefficient  combining  diffusion  and  hydrodynamic 
dispersion  (cm  /hour),  B is  the  average  pore  density  derived 
from  the  volumetric  water  content  and  C is  the  average  pore 
water  velocity,  and  D is  the  ratio  between  soil  bulk  density 
and  volumetric  water  content. 

In  a soil  system  where  the  terms  Bda/dx  and  Ddf/dx 
approach  zero  in  a soil  column  of  a semi-infinite  length, 
over  a small  time  t,  Eq.  (35)  can  have  the  following  solu- 
tion given  by  Lingstrom  et  al.  (1968): 


a(x,t)  = a erfc 


X 


2Kt 


(36) 


where  a^  is  the  initial  concentration,  and  K is  a constant 
depending  on  the  free  energy  of  adsorption. 
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Root  Effects 

The  effect  of  roots  in  a soil  solution  arises  because 
ions  movements  occur  into  and  out  of  the  roots,  with  the  net 
balance  being  an  influx.  This  uptake  is  mainly  associated 
with  ion  transport  by  diffusion  and  mass  flow  and  is  propor- 
tional to  the  concentration  at  the  root  surface.  It  is 
assumed  that  the  rate  of  plant  uptake  is  equal  to  the  rate 
of  loss  of  the  solute,  so  that  the  reaction  is: 


da 

dt 


= -6a 


(37) 


where  t is  the  reaction  time,  and  6 is  a constant.  A solu- 


tion of  the  above  differential  equation  is  given  by  Baldwin 
et  al.  (1973)  in  the  following  form: 


a 


aoexp 


-wt/(w,  + 


W2X 


In  ^) 


(38) 


X -r 

where  a^  is  the  initial  concentration  in  solution,  and  where 
w,  w^,  and  w^  are  constants,  r is  the  root  radiu,,  and  x is 
the  distance  from  the  root  surface.  This  equation  describes 
the  change  of  the  solution  concentration  at  a certain  dis- 
tance from  the  root  surface  with  time. 


CHAPTER  II 

KINETICS  OF  ADSORPTION  AND  DESORPTION 

The  interfacial  region  between  soil  colloid  and  solu- 
tion is  a center  of  intensive  chemical  and  physical  activi- 
ties. The  type  of  activities  retaining  our  attention  here 
is  the  adsorption  and  desorption  processes  of  ions. 

Adsorption 

The  adsorption  of  an  ion  on  the  colloid  surface  can  be 
considered  as  a second  order  reaction  involving  an  ion  (A) 
and  the  sites  (S)  on  the  surface: 

^1 

A + S i AS  (39) 

where  is  the  rate  constant  of  adsorption. 

By  assuming  that  the  reacting  sites  on  the  colloid 
surface  are  reacting  species  equivalent  to  (A),  then 

- ^ = k^[A]^  (40a) 

and  rearranging  this  is 

= /^k,  dt  (40b) 

[A]^  t=0  ^ 

or  can  be  written  as 

— — = k,t  (41) 

[A]  [A]^ 

The  plot  of  1/[A]  versus  time  t should  yield  a straight  line. 


26 


27 


Desorption 

The  desorption  of  any  adsorbed  ion  (AS)  can  be  con- 
sidered as  a pseudo-first  order  reaction: 


A + S 

where  ]<2  is  the 
tiation  of  Eq. 


rate  constant  of  desorption.  By  differen- 
(42),  the  following  can  be  obtained: 


k^CA]  (43) 

And  by  integrating  the  above  differential  equation,  this 
becomes 

[A]  = [A]  exp(k„t)  (44) 

o 2 

where  tA]^  is  the  initial  concentration  of  the  ion  specie  A, 
t is  the  reaction  time,  and  k2  is  the  rate  constant  of 
desorption . 


Adsorption  and  Desorption  Relationships 
Since  the  adsorption  and  desorption  phenomena  are 
simultaneous  a combination  of  the  Eq . (41)  and  Eq . (44)  is 
formed  as  follows : 

A + S ==P-==  AS  (45) 

^2 

d[A]/dt  = -k^[A]^  (46) 

d[A]/dt  = k2[AS] 

where  A,  S,  t,  k^ , and  k2  are  described  in  the  precedent 


(47) 
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sections.  By  dividing  the  Eq . (46)  and  Eq . (47),  and  rear- 
ranging, it  can  be  obtained: 

d[A]  _ ^ d[A] 
dCA]^  ^2  [A] 

1/[A]  - 1/[A]  = (k,/k„)ln([AS]/[AS]  (48) 

o 1 2 o 

If  we  set  a = 1/[A]  - 1/[A]  and  b = [AS],  then  Eq . (48)  can 

o ^ 

be  written  as 

^2 

b = b exp(:r—  a ) ( 49  ) 

o 

where  b is  the  amount  initially  adsorbed. 

The  above  equations  conbining  adsorption  and  desorption 
phenomena  are  valid  only  at  equilibrium,  there  is  need  to 
propose  another  model  combining  adsorption  and  desorption 
but  valid  at  any  reaction  time. 

A Kinetic  Model  for  Adsorption-Desorption 

Energy  Constant  Characteristic  of  the  Phase 

For  a collision  to  be  effective  in  producing  molecular 
species  from  reactants,  some  amount  of  energy  must  be 
available  to  allow  for  the  necessary  bonds  to  break  and  be 
formed.  In  1889,  Arrhenius  suggested  that  molecules  must 
get  into  an  activated  state  before  they  become  reactive.  In 
any  system,  an  equilibrium  exists  between  ordinary  and  active 
molecules  and  only  the  latter  are  rich  enough  energetically 
to  undergo  reaction 
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A ====  (A-B)^  ===  B (50) 

# 

where  (A-B)  is  the  spatial  configuration  of  the  transition 
state.  The  species  (A-B)  does  not  represent  the  active 
molecule  of  Arrhenius.  His  active  or  energized  species  are 

rather  a few  reactant  molecules  having  sufficient  energy  to 

• # 

get  into  the  transition  state,  (A-B)  , or  activated  complex 

condition,  but  not  necessarily  having  the  spatial  configura- 
tion corresponding  to  the  transition  state  (Eying  and  Eying 
1963). 

Let  us  consider  the  following  system 


/Timiin  / — 

////  C III  - B 

///I////////  - - 

7!  - / 

//I/////// 

1 

//I/////// 

1 

1 

1/1111111' 

/ 

///////// 

/ 

/ 

Figure  1.  Schematic  representation  of  solid-interface- 
solution  system. 

where  A is  the  solution  phase,  B is  the  interfacial  phase, 
and  C is  the  solid  phase.  Here  A,  B,  and  C are  different 
degrees  of  energized  species  of  the  same  ion  (or  molecule) 
with  respect  to  the  solid  surface;  these  are  considered 
characteristic  of  solution,  interfacial,  and  solid  phase, 
respectively.  The  surface  sites  at  which  ions  take  the 
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form  C are  at  random  on  the  solid  surface.  It  is  believed 
that  the  sites  associated  with  the  lowest  energy  are  the 
first  to  be  filled  by  the  ions  taking  the  form  C. 

Formulation  of  the  Model 

In  the  system  described  above,  it  can  be  assumed  that 
the  following  type  of  reaction  is  taking  place: 


A 


^-1 


k2 

B,  = - C 

b k_2  c 


(51) 


where  A,  B,  and  C are  as  described  previously  and  subscripts 
a,  b,  and  c designate  the  amounts  of  the  ion  in  the  forms  A, 
B,  and  C per  unit  volume  of  the  system.  The  terms  k^  k2  are 
the  rate  constants  for  forward  reaction  and  reverse  reac- 
tions, correspondingly  are  k-1  and  k-2. 

The  rate  of  a reaction  can  be  expressed  by  the  rate  of 
change  of  the  concentration  of  the  ion  species: 


da  , , , 

dt  - ^-1^  - ^1^ 


(52a) 


which  in  expanded  form  becomes 


dt  ^ ^1^  ■ k_ib  + k_2C  -k2b  = k^a  - (k_^+k2)b  + k_2C  (52b) 
and 


dc 

dt 


k2b  -k_2C 


(52c) 


Since  any  colloid  cannot,  within  a finite  period  of  time, 
indefinitely  fix  an  ion  in  the  C energy  form,  then  c must 
reach  a maximum  s at  equilibrium.  In  Eq . (52b),  it  is 
assumed  that  k_2  is  very  low  compared  to  k2. 
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The  simultaneous  solution  of  the  differential  for  Eq. 
(52a)  and  Eq.  (52b)  is 

w,t  w„t 

a = ae  and  b = Be  (53) 

where 


w^t  W^t  Wj^t 

aw^e  = -k^ae  + k_^Be 


(54a) 


and 

w,t  W-|t  w„t 

av]^e  = k^ae  - (k_^  + K2 ) Be  (54b) 

From  the  Eq.  (54a)  and  Eq.  (54b),  the  coefficients  a and  B 
can  be  determined  as  a function  of  the  rate  constants. 

In  any  soil- solution  system,  equilibrium  is  attained 
when  the  rate  of  change  of  each  ion  in  the  different  phase 
becomes  nil,  so  that  it  is  then 


da 

dt 


0 


such  that  k2b  - k_2C 
duced  as 


dc 

dt 


0 


0 and  c=s,  then  the  term  B is  de- 


Similarly , 
value,  a i 


„ ’'-2 

6 = 


if  k ^b  - k^a 


s found  to  be 


(55) 

0,  then  substituting  in  the  B 


a 


s 


e 


-w. 


t 

o 


In  the  initial  a and  b expressions,  shown  in  Eq . 


(SB) 
(53)  , 


it  can  be  shown  that 
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and 


k , k „ w, (t-t  ) 
^ ^ -1  -2  1 o 

^ ^ k — ic — ® 


k „ w^(t-t^) 

W — ^ ^ ^ 

b = s — e 
^2 


(57a) 


(57b) 


The  Eq.  (52c)  can  be  expressed  in  the  following  as  a 
first  order  linear  differential  equation: 


j w„t 

d|  " ^-2^  = e^2^ 


(58) 


The  solution  to  this  differential  equation  to  satisfy  the 
condition  where  c = 0 at  t = o can  be  written  in  the  form 


6k, 


c = 


k^+w^ 


”2''  k ,t 
e - e -2 


(59a) 


If  we  assume  that  k_2  is  negligible,  c can  be  further 
simplif ied 


6k 


c = 


k^+w^ 


2 "2'' 
e 


(59b) 


where  w^  and  w^  are  constants. 


CHAPTER  III 

THERMODYNAMICS  OF  ADSORPTION-DESORPTION  REACTION 


At  the  microregion  in  the  colloid  or  soil  solution 
extending  from  the  surface  to  the  outer  limit  of  the  first 
adsorbed  layer,  it  is  assumed  that  the  electrical  effect  due 
to  surface  charge  is  negligible  to  avoid  the  difficulty  in 
estimation  of  the  electric  field  effect.  In  this  microregion, 
there  is  ion  interaction  between  both  the  surface  and  the 
electric  potential  effects  (which  are  induced  by  the  surface 
charge ) . 

The  Surface  Charge 

The  colloidal  behavior  depends  on  how  the  surface 

charge  originated.  There  are  two  type  of  colloidal  charges. 

The  first  type  has  the  charges  due  to  crystalline  imperfec- 

4+  3 + 

txon,  such  as  isomorphic  substitution  of  Si  by  A1  or  by 

2 + 

other  cations  such  as  (Mg  ) having  a lower  charge.  This 
type  of  charge  is  found  in  the  crystal  lattice  of  clays  such 
as  montmorillonite  and  vermiculite.  In  these  cases,  the 
charge  density  is  constant  per  unit  of  surface  area.  The 
second  type  occurs  where  the  surface  charge  may  be  created 
by  preferential  adsorption  of  a certain  ion,  such  as  that  by 
hydroxyl  or  phosphate  ions.  In  this  case,  the  charge  arises 
at  the  exterior  edge  of  crystals  or  lattice,  thereby  inducing 
a constant  surface  potential  (van  Olphen  1977;  Bolt  1976). 
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Gouy-Chapman  and  Stern  Theories  of  the  Double-Layer 

From  Gouy-Chapman  theory,  the  charge  density  can  be 
deduced  considering  the  electroneutrality  condition  (van 
Olphen  1977;  Sennet  and  Olivier  1965): 


where  p is  the  space  charge  density  (or  net  sum  of  positive 
and  negative  ion  concentration,  o is  the  surface  charge,  and 
is  the  surface  potential.  The  Gouy-Chapman  theory  gives 
an  over-estimation  of  the  double  layer  capacity,  namely 


where  K is  the  reciprocal  of  the  double  layer  thickness,  a 
is  the  ion  concentration,  e is  the  electron  charge,  z is 
the  valence  of  the  ion,  and  e is  the  dielectric  constant. 

Stern  recognized  the  importance  of  the  ion  sizes  near 
the  surface.  He  proposed  that  the  counter  ions  could  be 
divided  between  a diffuse  layer  and  an  immobile  surface 
layer  of  thickness  6 able  to  contain  a maximum  number  of 
counter  ions  per  unit  of  surface  area.  This  may  be  ex- 
pressed as 


CO 

0 = -/^pdx 


(60) 


a 


K = (8z  e a/e  kT) 


(62) 


(63) 


where  o is  the  charge  corresponding  to  a monolayer  of 
m 
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coun'tGr’ions ) <f>  is  "the  van  der  Wall  Gnengy,  A is  ihG  frGquGncy 
factor,  1'6  is  the  GlGctric  potantial  at  thG  StGrn  Layar. 


Zgpo  Point  of  ChargG  (ZPC) 

On  a niGtal  oxida  surfaca,  chargas  ara  creatad  by  tha 
adsorption  and  dasorption  of  or  0H“  ions  which  are 
affected  by  their  concentrations  in  solution.  The  Parks  and 
da  Bruyn  (1962)  model  is 

(64) 


where  M is  A1  or  Fa  with  A and  C are  the  associate  anion 
and  cation,  respectively.  From  the  Eq . (64),  it  can  be  seen 
That  there  is  a pH  at  which  the  net  surface  charge  is  zeroj 
this  point  is  termed  the  zero  point  of  charge  (ZPC). 

For  metal  oxides,  such  as  descriWd  above,  H"*”  and  0H~ 
ions  are  potential  determining  ions . By  approximation  to 
the  Hernst  equation,  Keng  and  Uehara  (1974)  reported  that 

.+ 


RT  T H 

o 


(65) 


where  R is  the  gas  constant,  T is  the  absolute  temperature,  F 
is  the  Faraday  constant,  and  <<25  mV;  the  Gouy-Chapman 
double  layer  equation  is  reduced  according  to 
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% = = W (0-059)(pH^  - pH)  (66) 

Thermodynamics  of  Adsorption 

Effect  of  the  Electric  Field 

In  this  system  introduced  in  Fig.  (1),  the  compartments 
A,  B,  and  C can  be  considered  as  being  three  phases  of  a 
particular  ion  at  different  potential  energy  levels.  In 
each  phase,  molecules  are  at  different  energy  levels  and 
varying  within  a range  of  potential  characteristics  of  each 
phase . 

Around  a colloid  center,  the  ions  exhibit  a Boltzman 
type  distribution: 

c = b exp  [-F(T^  - 'F^)/RT]  (67a) 

b = a exp  C-F(Tj^  - T^)/RT]  (67b) 

Supposing  that  1'  = 0,  it  can  be  shown  that 

a 

''o  = '•'o  = - r I 

where  T is  the  surface  potential,  and  T are  electric 
o ^ b a 

potentials  that  is  characteristic  of  the  phases  B and  A, 
respectively . 

Free  Energy  as  a Function  of  Distance 
From  the  outer  limit  of  the  Stern  layer,  both  Gaussian 
and  Gouy-Chapman  concepts  can  be  applied  to  a particular 
ion.  The  movement  of  an  ion  (or  molecule)  results  from  a 
succession  of  collisions  that  may  move  it  at  random  in  a 
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positive  or  negative  direction.  These  ionic  concentrations 
follow  the  Gaussian  distribution: 

a = 7 » exp(-x^/4Dt)  (68) 

2(nDt)  ' 


where  M is  the  amount  of  substance  deposited  at  the  plane 
when  X = 0 at  the  time  t = 0,  and  D is  a constant.  Assuming 
that  the  Gouy-Chapman  distribution  holds  in  the  following 
equation : 

T = exp(-Kx)  (69) 

where  the  free  energy  can  be  expressed  as 

G.  ^ = G°  + zFT  + RTln  a (70a) 

and 

* zFT^exp(-Kx)  + RTln  — 


At  constant  time  t,  this  expression  becomes 

RT 


dG  ,,  -Kx 

T—  = -zFT  Ke 
dx  o 


2W 


X 


(71) 


The  limit  between  the  physically  adsorbed  layer  and  the  com- 
pletely free  ion  in  the  phase  A may  be  defined  when  dG/dx=0. 
If  we  set  y = -Kx  and  y = -RT/ 2Dt/zFT^K,  then  the  expression 
can  be  written 


y = 0 


(72) 


and  become 


y^  + 2y(l  - h 


+ 1 = 0 


where 
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(e^  = 1 + y + 1/2  + . . . ) 

V7hen  this  quadratic  equation  gives  two  roots , they  may  be 
indicators  of  the  transition  states  between  different 
degrees  of  adsorption. 

Free  Energy  for  Irreversible  Fixation 
The  active  transport  of  an  ionic  substance  against  the 
gradient  potential  is  determined  by  the  difference  between 
the  total  potential  within  each  phase  (A,  B,  and  C).  The 
maximum  energy,  other  than  expansion  work  for  a change  of 
ion  activities,  is  expressed  by  at  constant  temperature 

and  pressure.  The  transfer  of  ions  between  the  phase  B and 

C is  due  to  a potential  energy  minimum  for  phase  C written 

c B 

as  G.  and  maximum  for  phase  B written  as  G.  . In  this 
1 ,m  ^ 1 ,m 

case 

g9  = z.F4'^  + X.  (73) 

1 ,m  1 1 

and,  correspondingly, 

G?  = z.FH'®  + RTln  b (74) 

1 ,m  1 

where  G^  ^ is  also  the  Gibbs  free  energy  of  the  ion  i per 
unit  mole,  is  defined  as  the  minimum  energy  for  irrevers- 
ible fixation  (including  only  the  electro-chemical  energy), 

B C 

F is  the  Faraday  constant,  T and  T are  the  electric  poten- 
tial in  the  phases  B and  C,  and  the  term  RTln  b is  the  chem- 
ical potential.  The  transfer  of  ions  between  the  phases  B 
and  C is  governed  by  the  difference  in  total  energy  between 
the  two  phases. 
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k_2 

AG.  = -46az.F  + RTln  s ^ + RTw„(t  - t ) - X.  (75) 

1 ^2  z o 1 

where  'i'  = = 4II6a  which  is  the  Gaussian 

equation  for  a molecular  condenser.  This  by  substitution 

gives 

2 

AG.  ^ = 4n6az.F  + RTln  s t— ^ - X.  (76) 

i,m  1 )<2  1 

When  the  conditions  are  at  equilibrium,  or  AG  = 0 and  t = t^, 
then 

k 

X.  = - 4n6az.F  + FTln  s (77) 

1 1 K2 

It  is  assumed  that  the  minimum  energy  X^  is  independent 

of  the  state  of  equilibrium,  then  for  a given  colloid  surface 

which  has  a particular  surface  charge,  the  free  energy  change 

AG.  is  reduced  to 
1 , m 

AG.  ^ = RTw„(t  - t^)  (78) 

1 , m / o 


Relationship  of  Surface  Charge  to  Surface  Potential 
During  the  overall  transfer  of  ions  from  the  phase  A 
to  phase  C at  equilibrium  conditions,  after  substituting  the 
X^  value  of  Eq.  (77)  and  assuming  H'  =0,  the  relationship 
becomes 


d>F 


o 


446 
a . F 


■da 


1 

where  6 is  the  Stern  layer  thickness, 
F is  the  Faraday  constant,  and  a^.  is 


o 

the 


(79) 

is  the  surface  charge, 
ion  valence. 
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Surface  Tension  and  Specific  Adsorption 
In  the  interfacial  colloid-liquid  region,  some  ions 
may  enter  into  the  coordination  shell  of  surface  atoms.  As 
a result,  there  is  a modification  of  the  colloid  surface  (A) 
(expansion  or,  contraction) . This  work  required  to  expand 
or  contract  the  surface  when  divided  by  the  change  in  sur- 


face  area  is  the  surface  tension  (y)  (Atkins  1978). 
following  series  demonstrates  the  relationship: 

The 

g9  = z.FTc  + Ay 

(80a) 

g9  = z.FTc  + X.c 
11  1 

Differentiating  the  Eq.  (80a)  and  Eq . (80b)  gives 

(80b) 

dG^*“  = z^FTdc  + z^FcdT  + ydA  + Ady 

(81a) 

dG9  = z.FTdc  + z.FcdT  + X . dc  + cdX . 
11  1 1 1 

However,  it  is  known  that 

(81b) 

dG9  = z.F'Fdc  + X.dc 
11  1 

(82a) 

dG9  = z.FH'dc  + ydA 
1 1 ' 

(82b) 

By  comparing  Eq . (81a)  to  Eq.  (82b)  and  Eq . (81b)  to  Eq . 


(82a)  it  follows  that 
cdXE  = Ady 

(83) 

The  differentiation  of  Eq . (77)  gives  dX^  = -4n(Sz^F.do,  so 


it  can  be  deduced  that 

, -4n6z. 

dy  _ 1 _ 

do  A 

(84) 

CHAPTER  IV 

MATERIALS  AND  METHODS 
Goethite  Preparation 

Three  goethite  preparations  were  made  by  mixing  sepa- 
rately an  equal  volume  (50  ml)  of  1 N FeCCD^  with  a similar 
volume  of  2 N NaOH,  1 N NaOH,  or  0.5  N NaOH  in  order  to  have 
OH/Fe  mole  ratios  of  6,  3,  and  1.5,  respectively,  in  the 
mixtures.  The  suspensions  were  allowed  to  age  for  one  week 
to  induce  crystal  formation  and  growth.  The  addition  of 
phosphorus  to  the  goethite  preparation  was  made  for  P/Fe  = 3 
ratio  both  at  the  beginning  and  by  adjustment  at  the  end  of 
ageing  process  (one  day  of  phosphorus  reaction  is  allowed). 
In  order  to  examine  the  effects  of  phosphorus  concentrations 
on  goethite  crystallization,  different  levels  of  phosphorus 
in  2 N NaOH  were  mixed  with  1 N Fe(Cl)2  a.nd  allowed  to  age. 
After  one  week  of  ageing  the  suspension  at  either  room 
temperature  or  at  55°C,  the  samples  were  washed  by  dialysis 
against  distilled  water  during  one  week;  the  water  was 
changed  after  intervals  not  greater  than  12  hours.  At  the 
end  of  the  ageing  process  and  washing,  the  samples  were 
freeze-dried . 

The  existence  of  hydroxyl  def ormational  vibrations  in 
the  region  1200  to  1000  cm“^  was  investigated  by  two  methods 

i 
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1)  goethite  and  phosphated  goethite  were  digested  in  D2O  for 
24  hours  to  replace  the  surface  OH  by  OD , or  2)  phosphated 
goethite  samples  were  equilibrated  with  different  salts 
(0.1  N KCl,  0.1  N KNO^,  and  0.1  N NajSO^)  and  water  to 
displace  surface  phosphates.  The  samples  were  separated 
from  solution  and  dried  for  the  infrared  spectroscopy 
studies . 

An  attempt  was  made  to  prepare  FeOOD  by  digesting  three 
times  0.9  g of  FeC1^.6Fl2)  in  D2O  and  drying  the  suspension 
in  order  to  make  FeCl2.6D20.  The  residues  were  dissolved  in 
10  ml  of  D2O  and  mixed  with  10  ml  2 M NaOD;  the  suspensions 
were  then  aged  and  dried  as  described  above. 

Infrared  and  X-Ray  Diffraction  Techniques 

One  milligram  of  the  freeze-dried  sample  was  mixed  with 
400  mg  KBr  to  make  pellets  samples  used  for  the  infrared 
study.  The  infrared  spectra  are  obtained  by  using  Perkin- 
Elmer  567  grading  infrared  spectrophotometer. 

Prior  to  X-ray  diffraction  analysis  using  a general 
electric  XRD700  instrument,  and  thin  films  of  goethite 
samples  were  made  on  glass  slides  and  allowed  to  air  dry. 

Some  Factors  Affecting  the  Kinetics  of  Adsorption-Desorption 
The  adsorption  or  desorption  isotherm  is  obtained  by 
plotting  the  amount  of  ion  adsorbed  or  desorbed  per  unit 
v;eight  of  adsorbent  versus  the  solution  concentration.  The 
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solid  and  the  solution  were  equilibrated  through  agitation 
for  selected  periods  of  time  and  temperatures.  The  equilib- 
rium solution  was  removed  after  centrifugation  and  the 
phosphorus  concentration  determined  through  blue  color 
development  as  ascorbic  acid  molybdophosphoric  complex. 

Effects  of  pH 

The  adsorption  for  different  periods  of  times  (0  to  72 
hours)  was  done  by  equilibrating  goethite  or  soils  with 
phosphorus  solution  (Syg  P/ml)  or  (20  yg  P/ml)  adjusted  to 
pH  2,  4,  8,  9.5,  and  11.  Solid  phase  was  separated  by 
centrifuging  at  5000  revolutions  per  minute  (rpm) . Phosphate 
desorption  was  carried  out  on  samples  which  reached  the 
equilibrium  state.  The  desorption  was  conducted  in  water, 
adjusted  to  the  above  pH,  during  20  minutes,  6 and  12  hours. 
The  pH  of  the  solution  was  measured  by  placing  the  glass 
electrode  in  the  supernatant  after  centrifugation. 

Effects  of  the  Supporting  Electrolytes 

The  effects  of  the  type  of  cations  on  the  adsorption 
were  studied  by  mixing  20  mg  of  goethite  with  20  ml  each  of 
the  phosphate  solutions  containing  1,  3,  5,  and  10  yg  P/ml 
in  0.01  M salts  of  NaCl,  CaCl2,  and  AlCl^,  each  in  separate 
experiments.  With  soils,  the  phosphate  solution  concentra- 
tions used  were  0,  10,  20,  40,  and  80  yg  P/ml.  The  extent 
that  each  of  the  cations  (Ca,  Ha,  and  A1 ) contributed  to  the 
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degree  of  reversibility  of  the  adsorbed  phosphate  was  studied 
by  conducting  the  desorption  with  solution  of  the  correspond- 
ing electrolyte  (as  in  the  adsorption)  but  without  phosphate. 
A comparison  of  the  phosphate  desorption  by  different  anions 
in  various  salts  (KCl,  K^SO^,  KNO^,  KCIO^,  and  NaOH)  was 
made  for  goethite-P  solution  (one  gram  of  solid  with  one 
liter  of  salt  solution  where  the  goethite  had  previously 
been  treated  with  5 ug  P/ml). 

Time  Effects  on  Adsorption  and  Desorption 

Adsorption  of  phosphates  on  goethite-solution  containing 
0,  1,  3,  5,  7,  and  10  yg  P/ml,  and  on  soil-solution  contain- 
ing 0,  5,  10,  20,  40,  and  80  yg  P/ml  was  investigated  after 
equilibration  during  different  periods  of  time  ranging  from 
1/10  hour  to  80  hours.  At  the  end  of  each  reaction  time, 
the  solid-solution  was  centrifuged  and  the  phosphate  concen- 
tration remaining  in  the  supernatant  solution  determined. 

The  solid/solution  ratio  was  1/1000  for  goethite-solution 
and  1/20  for  soil-solution.  Kinetic  studies  of  phosphate 
adsorption  were  done  at  different  conditions:  1)  when 

different  types  of  electrolytes  were  used,  2)  for  one  elec- 
trolyte (CaCl2)  at  different  concentrations  (.01  N Ca, 

0.1  N Ca,  1.0  K Ca),  and  3)  at  pH  2 and  pH  10.  After  the 
adsorption  proceeded  until  the  equilibrium  state  v;as  reached, 
desorption  was  conducted  during  20  minutes,  6 and  12  hours 
in  order  to  determine  the  minimum  time  required  for  complete 
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desorption  of  this  form  of  phosphate.  This  minimum  time  was 
also  employed  during  desorption  of  phosphate  when  using 
either  a desorbing  solution  at  selected  pH  (2  to  11)  or  for 
different  electrolytes.  For  those  soils  where  it  was  observed 
that  phosphate  sorption  did  not  follow  the  Langmuir  isotherm, 
an  attempt  was  made  to  determine  the  extent  of  changes  in 
aluminum  and  iron  phosphate  by  using  the  method  of  Peterson 
and  Corey  (1966).  For  aluminum  phosphate,  one  gram  of  soil 
sample  was  washed  with  2 N NaCl  to  remove  exchangeable 
cations,  the  suspension  was  centrifuged  and  the  supernatant 
solution  discarded.  Then  20  ml  of  0.5  N NH^^F  at  pH8.2  were 
added  and  the  suspension  was  shaken  for  an  hour  and  then 
centrifuged  for  the  P determination.  For  the  iron  phosphate 
fraction,  the  samples  were  washed  with  2 N NaOH  solution, 
centrifuged,  and  the  supernatant  retained  for  P 
determination . 

Surface  Charge  as  Affected  by  Phosphate  Adsorption 

The  method  of  Lavardiere  and  Weaver  (1977)  was  used  to 
determine  the  net  electric  charge . The  procedure  was  to  add 
20  ml  each  of  0.01  N CaCl2  and  1.0  N CaCl2  solutions  per  one 
gram  of  soil  sample  or  20  mg  of  goethite.  Subsequent  titra- 
tions of  the  suspensions  were  made  with  0.01  N HCl  or  0.01 
N NaOH,  by  adding  0.1  to  0.3  ml  at  a time  from  a microburet 
at  2-minute  interv^als.  Continuous  stirring  was  maintained 
and  the  pH  read  before  the  addition  of  either  base  or  acid. 
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A blank  titration  was  made  for  the  same  volume  of  CaCl2 
solution.  The  amount  of  H or  OH  adsorbed  at  a given  pH 
was  calculated  by  the  difference  between  the  amount  of  H^  or 
OH  added  and  that  required  to  bring  the  blank  solution  of 
the  same  volume  and  salt  concentration  to  the  same  pH  as 
the  soil  or  goethite. 

Soil  pH,  Iron  Oxide,  and  Extractable  P 

The  soils  used  were  obtained  from  Georgia  (Cecil,  Ap 
horizon),  Colorado  (roadside  cut  near  Ft.  Collins),  and  Kenya 
(latosol,  sampled  at  15-30  cm,  near  Kabete). 

The  pH  was  measured  in  1:1  soilto  solution  suspension 
for  10  g of  soil  with  10  ml  of  H^O  or  10  ml  of  1 N HCl , using 
a glass  electrode.  Iron  oxide  was  determined  by  the 
dithionite-bicarbonate  extraction  and  colorimetric 
determination  of  Fe  as  the  ferrous  orthophenonthroline  complex. 
Extractable  P was  determined  by  three  of  the  methods  outlined 
by  Ballard  (1979).  These  data  are  shown  in  Table  1. 
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Table  1. 

Soil  pH,  iron  oxide 

and  extractable 

P. 

Soils 

Georgia 

Colorado 

Kenya 

pHCH^O) 

5 . 88 

8.17 

6.20 

pH(KCl) 

4 .95 

7 .44 

5 .75 

% Fe202 

2 . 24 

0.67 

3.95 

P,  ppm 

0 . 05 

N HCl  +0.025  N 

5.60 

0 . 30 

0.40 

P,  ppm 

0.03 

N NHj^F  + 0.1_N 

HCl* 

13.00 

7.40 

1.40 

P,  ppm 

0.5  M NaHCOg-' 

2 .80 

1.20 

1.60 

Reagents  described  by  Ballard  (1974) 


CHAPTER  V 

RESULTS  AND  DISCUSSION 

Goethite  and  Phosphated  Goethite  Studies 
by  InfraRed  Spectroscopy 

After  the  mixing  of  FeCI^  and  NaOH  solutions,  some 
precipitates  formed,  but  the  particles  formed  at  precipita- 
tion were  not  yet  crystalline.  The  changes  in  precipitates 
to  crystals  originated  from  a discrete  center  or  crystal 
nuclei  (twinned  and  acicular  crystals)  produced  by  different 
mechanisms  (Atkinson  et  al.  1968).  The  conditions  governing 
the  formation  and  nature  of  these  crystal  nuclei  are  not 
well  known.  Apparently,  the  number  of  crystal  nuclei  was 
increased  by  both  temperature  and  hydroxyl  ion  concentrations. 
Iron  and  hydroxyl  ions  would  be  attracted  to  the  centers  of 
crystal  growth  as  they  lose  their  energies.  As  iron  and 
hydroxyl  ions  were  involved  in  the  formation  of  the  crystal, 
ageing  reactions  follow  the  sequence  of  lepidocrocite  to 
goethite  (Murphy  et  al.  1975).  Another  way  would  be  that 
ferrhydrite  forms  first  or  just  goethite  was  formed.  The 
growth  of  the  crystal  occurred  as  the  ions  change  phase 
during  deposition  from  the  solution.  Hsu  (1972)  found  that 
the  removal  of  hydroxyl  ions  from  the  solution  resulted  in  a 
drop  in  pH  value  of  the  solution  as  the  time  increased. 
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Since  the  ions  are  changing  phase  from  the  solution  (higher 
degree  of  randomness)  to  the  solid  phase  (lower  degree  of 
randomness),  the  entropy  change  during  the  process  from 
solution  to  crystalline  phase  must  be  negative.  The  entropy 
of  ions  within  the  crystal  would  be  lower  than  those  at  the 
crystal  surface.  From  infrared  spectra,  Kiselev  and  Lygin 
(1975)  believed  it  was  possible  to  estimate  from  adsorption 
entropy  what  was  the  degree  of  freedom  of  the  adsorbed 
molecules . 

Goethite  Structure  Identification  by  Infrared 

The  4000  - 2000  cm  ^ Region 

The  goethite  structure  was  interpreted  through  the 
identification  of  OH  and  FeO  vibrations.  The  hydroxyl  ion 
vibrations  occurred  at  two  strong  bands  in  the 
3700  - 2000  cm  ^ region,  centered  at  3400  cm~^  and  3200  cm~^ 
(see  Fig.  2 and  Fig.  3).  According  to  Nakamoto  (1978), 
lattice  water  absorbed  at  3550  - 3200  cm~^,  he  reported 
that  the  strong  band  centered  at  3400  cm~^  was  due  to  anti- 
symmetric and  symmetric  OH-stretching  of  water.  It  could 
easily  be  recognized  that  the  3400  cm~^  band  was  not  charac- 
teristic of  goethite  crystal  structure  since  this  band 
appears  even  when  no  goethite  is  present,  as  determined  by 
X-ray  diffraction  which  showed  no  crystalline  material.  It 
is  of  interest  to  observe  that  the  intensity  of  the  3200  cm"^ 
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3500  3000  2500 


Frequency  (cm 

Fig.  2.  Infrared  bands  of  (A)  goethite ; (B)  phosphate  added 
to  goethite  at  the  end  of  ageing;  and  (C)  phosphate 
added  to  goethite  at  the  beginning  of  ageing. 

Curves  A,  B and  C are  for  22°  C and  curves  A',  B', 
and  C are  for  55°  C. 
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Infrared  bands  of  (A)  goethite ; (B)  phosphate  added 
to  goethite  at  end  of  ageing;  (C)  phosphate  added 
to  goethite  at  beginning  ageing.  Curves  A,  B,  and 
C are  for  22°  C,  and  curves  A’,  B',  and  C are  for 
55°C. 


Fig.  3. 
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band  decreased  as  the  OH/Fe  ratio  in  the  ageing  solution  de- 
creased. The  band  centered  at  3200  cm”^  appeared  only  if 
goethite  structure  was  present.  As  a result,  it  was  con- 
cluded that  the  3200  cm  ^ band  corresponded  to  structural 
Fe  - OH  stretching. 

The  2000  - 3000  cm  ^ Region 

After  goethite  was  freeze-dried,  not  all  the  adsorbed 
water  was  removed.  Nakamoto  (1978)  indicated  that  the 
relative  velocity  of  the  oxygen  nucleus  compared  to  that  of 
hydrogen  nucleus  is  small.  This  meant  that  the  surface 
binding  of  water  through  the  oxygen  atom  to  goethite  would 
not  induce  a significant  change  in  the  overall  water  vibra- 
Tion  which  could  be  observed  at  around  1620  cm  ^ as  depicted 
for  water  structure  below 


/k 

H H 

d N 


/\ 

H H 


H H 


Case  V. 


Case  V, 


Case  V. 


The  above  three  normal  modes  of  vibration  in  H^O  are  infra- 
red active.  The  bending  vibration  ^ is  centered  at  1620 
cm  ^ and  the  water  stretching  bands  (v^  and  v^)  vibrate  in 
the  3400  cm 


region . 
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The  1300  - 700  cm  ^ Region 

The  two  main  bands  at  890  and  790  cm"^  would  be  assigned 
to  Fe-OH  bending  vibration  of  the  structural  hydroxyls. 

Busca  et  al.  (1978)  supported  this  finding  by  observing  that 
the  structural  OH  in  a-FeOOH  disappeared  upon  heating  to 
form  a-Fe^Og.  The  appearance  of  the  890  and  790  cm~^  bands 
always  indicated  goethite  crystallization,  which  can  be 
observed  in  Fig.  3,  4,  5,  6,  7,  9,  10  and  11.  From  the 
above  mentioned  figures,  it  was  observed  that  for  each  OH/Fe 
ratio,  the  way  in  which  phosphate  was  added  had  an  effect  on 
the  degree  of  goethite  crystallization.  From  Fig.  5,  as 
reported  by  Farmer  and  Palmieri  (1975),  typical  goethite 
could  be  identified  by  the  Fe  - OH  bending  vibration  at 
890  cm~^  and  790  cm"^.  In  Fig.  3,  where  OH/Fe  = 6,  there 
are  also  strong  bands  at  890  and  790  cm~^  in  all  cases, 
except  when  phosphate  is  added  at  the  beginning  of  ageing  of 
the  suspension  at  55°C.  In  the  latter  case,  the  increase  of 
temperature  to  55°C  favored  preferential  phosphate  binding 
to  the  iron  which  was  observed  by  the  strong  vibration  at 
1000  cm  As  shown  in  Fig.  3c,  phosphate  was  within  the 
goethite  structure  because  there  was  both  the  P - 0(Fe) 
vibration  at  1000  cm  ^ and  the  surface  binuclear  (FeO^POOH 
vibrations  at  1190,  1100,  and  1030  cm“^.  When  CH/Fe  = 3,  it 
was  not  possible  to  have  phosphate  v;ithin  the  crystal  because 
there  was  no  crystallization  at  room  temperature.  However 


TRANSHITTAKCE 


54 


Frequency  (cm 

Fig.  4.  Infrared  bands  of  (A)  goethite ; (B)  phosphate  added 
to  goethite  at  end  ageing;  and  (C)  phosphate  added 
to  goethite  at  beginning  ageing.  Curves  A,  B,  and 
C are  suspensions  at  22°C,  and  curves  A',  B',  and 
C are  for  suspensions  at  55°C. 
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Fig.  5. 


After_ Farmer  and_Palmieri  (1975).  Infrared  bands  of 
Goethite  and  lepidocrocite . 


TRANSMITTANCE 


56 


Fig.  6.  Infrared  bands  of  (A)  goethite;  (B)  phosphate  added 
to  goethite  at  end  ageing;  and  (C)  phosphate  added 
to  goethite  at  beginning  ageing.  Curves  A,  B and 
C are  for  suspensions  at  22°C  and  curves  A',  B', 
and  C are  for  55°C. 
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Fig.  7,.  Infrared  bands  of  (A)  Fe  hydroxide  material;  (3) 

phosphate  added  to  Fe  hydroxide  material  at  end  of 
ageing;  and  (C)  phosphate  added' to  Fe  hydroxide 
material  at  beginning  ageing.  Curves  A,  B,  C are 
for  suspensions  at  22°  C and  curves  A’,  B’ , C 
are  for  suspensions  at  55°  C. 
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there  was  weak  crystallization  at  55°C  that  favored  some 
phosphate  surface  binding  (see  Fig.  6b).  In  OH/Fe  = 1.5 
suspension,  after  one  week  of  ageing  either  at  room  tempera- 
ture or  at  55°C,  no  goethite  formation  was  observed.  There 
was  a single,  strong  vibration  at  1030  cm  ^ when  phosphate 
was  added  at  the  end  of  the  ageing,  but  the  vibration  was  at 
1000  cm  ^ when  the  phosphate  was  added  at  the  beginning  of 
ageing.  The  1030  cm~^  vibration,  which  is  the  P - OH  bending 
vibration,  indicated  the  excistence  of  phosphate  at  the 
surface,  and  the  presence  of  the  1000  cm~^  vibration  from 
the  P - 0(Fe)  showed  phosphate  directly  bound  to  iron  when 
phosphate  was  added  just  prior  to  the  ageing  process  (Parfitt 
et  al.  1975). 

Goethite  Identification  by  X-Ray 

The  X-ray  diffraction  peaks  for  goethite  were  at  4.19  8, 
2.70  8,  and  2.45  In  all  cases,  the  4.19  X is  the  most 

intense  peak  while  those  at  the  2.70  ^ and  2.45  K were  weak. 
Fig.  8.  The  presence  of  the  2.70  8 spacing  indicated  that 
some  hematite  might  be  present.  However,  because  the  2.45  8 
peak  was  present,  it  was  believed  that  considerable  amounts 
of  goethite  existed  (Schwertman  and  Taylor  1977).  According 
to  their  studies,  it  would  be  possible  for  hematite  to  exist 
along  with  goethite  because  the  hydration  of  hematite  would 
yield  goethite  with  a standard  free  energy  (AG°)  of  the 
reaction  varying  from  -0.2  to  0.4  kcal/mole. 


Angstrom  Spacing  (8) 
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Fig.  8.  X-ray  diffraction  patterns  of  (A)  goethite,  (B)  phosphated 
end  of  ageing,  and  (C)  phosphated  goethite  at  beginning  of 
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Factors  Affecting  Goethite  Crystallization 
Effect  of  OH/Fe  Ratio 

The  above  results  indicated  that  the  bands  at  3400  cm” 
3200  cm  890  cm  and  790  cm  ^ were  characteristic  of 
the  appearance  of  goethite.  For  the  same  ageing  period,  an 
increase  in  OH/Fe  favored  goethite  crystallization.  When 
the  OH/Fe  is  6,  the  bands  at  3200  cm”^,  890  cm”^,  and  790  cm 
were  strong,  indicating  that  goethite  structure  was  well- 
formed.  Where  OH/Fe  is  3,  after  a week  of  ageing,  goethite 
structure  was  apparently  present  only  if  the  suspension  was 
kept  at  55°C,  even  then  the  degree  of  crystallization  was 
less  intense  at  OH/Fe  = 3 than  it  was  where  OH/Fe  = 6.  This 
weak  crystallization  was  suggested  by  the  weak  band  at 
3200  cm  ^ (which  took  the  form  of  a shoulder)  and  by  the 
less  pronounced  intensities  of  the  890  cm  ^ and  790  cm  ^ 
bands.  However,  when  OH/Fe  is  1.5,  there  was  no  goethite 
crystallization  even  when  the  suspension  was  aged  at  55°C. 
The  above  conclusion  disagreed  with  that  of  Atkinson  et  al. 
(1974)  who  claimed  made  goethite  was  in  a suspension  at  28°C 
where  OH/Fe  is  2.0  when  aged  for  50  hours. 

Fhosphated  goethite 

The  application  of  phosphorus  weakened  the  goethite 
structure  because  some  phosphate  was  probably  incorporated 
within  the  lattice.  In  Fig.  4 and  Fig.  6 where  OH/Fe  is  3, 
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the  addition  of  phosphate  to  goethite  suspension  resulted  in 
bands  at  3200  cm  890  cm  ^ and  790  cm  ^ which  were  weaker 
than  those  for  OH/Fe  = 6 after  ageing.  When  OH/Fe  = 6 in 
the  suspension,  hydroxyl  ion  concentration  was  high  enough 
to  form  the  necessary  bonds  for  the  goethite  structure  to 
appear.  The  presence  of  phosphate  on  the  goethite  surface 
could  be  recognized  by  the  appearance  of  bands  in  the  1200  - 
1000  cm  ^ region,  but  vibrational  hydroxyl  deformations  also 
could  occur  in  this  same  region.  Parfitt  (1979)  stated  that 
the  P = 0 bond  had  stretching  vibrations  in  the  1190  cm  ^ and 
the  1030  cm  ^ region. 

When  the  phosphate  was  added  at  the  end  of  the  goethite 
ageing,  there  was  an  appearance  of  P = 0 vibration  at 
1190  cm  ^ and  the  1030  cm  ^ due  to  P - OH  vibration  in 
agreement  with  work  by  Parfitt  (1979).  At  high  temperature 
(55°)  even  when  OH/Fe  = 6,  the  phosphate  is  preferentially 
bonded  to  the  iron,  thereby  reducing  the  capacity  of  hydrox- 
yls ions  to  bind  freely  with  iron  to  form  goethite.  As  a 
result,  there  was  direct  binding  of  phosphate  to  iron  which 
was  confirmed  by  the  strong  band  at  1000  cm~^  assigned  to 
P - 0(Fe)  and  by  the  lack  of  the  Fe  - OH  bending  vibrations 
at  890  cm”^  and  790  cm~^.  This  is  confirmed  by  spectra 
given  in  Fig.  3.  In  Fig  9,  Fig  10  and  Fig  11  where  suspen- 
sion the  OH/Fe  ratio  was  either  6,  3,  or  1.5,  the  increase  in 
phosphate  concentrations  relative  to  the  iron  (P/Fe  = 0.032, 
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Fig.  9.  Infrared  bands  of  phosphated  goethite  at  beginning 
of  ageing  for  suspensions  of  various  P/Fe  values  at 
OH/Fe  = 6. 
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. Infrared  bands  of  phosphated  goethite  at  beginning 
of  ageing  for  various  P/Fe  ratios  at  OH/Fe  = 3.0. 
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Fig.  11.  Infrared  bands  of  phosphated  goethite  at  beginning 
of  ageing  for  various  F/Fe  ratios  at  OH/Fe  = 1.5. 
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0.32,  3.2)  induced  an  increase  in  the  intensity  of  P - 0(Fe) 
and  P - OH  stretching  vibration  at  1000  cm~^  and  1030  cm~^, 
respectively.  When  P/Fe  was  greater  or  equal  to  0.3,  these 
two  bands  overlapped,  showing  only  one  very  strong  band  at 
1000  cm~^. 

In  order  to  identify  the  vibration  of  hydroxyl  deforma- 
tion assumed  to  be  in  the  1200  - 1000  cm"^  region,  two 
methods  were  used:  1)  goethite  and  phosphated  goethite  were 

digested  in  D2O  for  24  hours  and  2)  surface  phosphate  was 
desorbed  by  different  anions  and  water.  After  the  samples 
were  dried,  the  infrared  spectra  were  similar  to  those 
shown  in  Fig.  12  and  Fig.  13.  The  weak  bands  due  to  OH 
deformation  were  displaced  by  OD  and  only  the  P - 0(Fe) 
stretching  vibration  at  1000  cm~^  and  those  for  the  Fe  - OH 
bending  at  890  cm  ^ and  790  cm  ^ persisted.  Evidently,  it 
could  be  assumed  that  the  bands  at  890  cm”^  and  790  cm~^ 
were  due  to  structural  Fe  - OH  bending  which  cannot  be 
affected  by  digestion  in  as  long  as  the  initial  goethite 

maintained  its  structure.  In  a further  study,  the  evidence 
for  hydroxyl  deformations  was  examined  when  the  phosphated 
goethites  were  desorbed  by  different  salts  (KCl,  KNO^, 
Na2S0^^)  and  water.  After  either  0.1  N KCl  solution  or  water 
desorption,  hydroxyl  deformations  were  not  removed  in  the 
1200  cm  ^ - 1000  cm  ^ region.  With  0.1  N Na2S0^  solution 
used  for  desorption,  the  infrared  spectra  showed  strong 
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Fig.  12.  Infrared  bands  of  1)  goethite  digested  in  D„0  and 
2)  phosphated  goethite  digested  in  D^O. 
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Fig.  13.  Infrared  bands  of  phosphated  goethite  after  desorp- 
tion by  (1)  0.1  N KCl,  (2)  H„0,  (3)  0.1  N KNO,, 

(4  ) 0.1  N Na^SO^^. 
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2-  -1 

bonds  of  SO^  ions  in  the  1200  - 1000  cm  region.  When 

0.1  N KNOg  was  used  for  desorption,  NO”  appeared  to  displace 

some  surface  phosphates  so  that  the  hydroxyl  vibration  in 

that  region  was  reduced  (Fig.  13).  It  can  be  said  that  NO^ 
2-  . 

and  ions  have  the  ability  to  desorp  readily  displace- 

able phosphate  at  the  surface.  However,  sulfate  ions  could 
not  be  used  to  provide  evidence  for  this  type  of  phosphate 
desorption  by  infrared  since  sulfate  ions  have  strong 
vibration  in  the  same  region  as  that  of  surface-bound  phos- 
phate. Parfitt  et  al . (1975)  assigned  the  appearance  of 
weak  bands  in  the  1200  - 1000  cm  ^ region  for  non-phosphated 
goethite  to  the  Fe  - OH  deformat ional  vibration  which  can 
partially  overlap  the  region  of  P = 0 stretching  and  P - OH 
bending  vibration. 

The  question  arose  about  the  type  of  hydroxyl  groups 
that  are  displaced  when  phosphates  are  added  to  the  goethite 
suspension.  Russel  et  al . (1974)  indicated  that  there  were 

three  types  of  hydroxyl  groups:  1)  where  OH  is  singly 

coordinated  to  an  Fe  atom  with  hydrogen  bond  interaction 
with  another  atom,  2 ) where  OH  is  coordinated  to  tv/o  Fe 
atoms  and  3)  when  OH  is  singly  coordinated  to  an'  Fe  atom. 
Parfitt  et  al . (1975)  considered  that  only  type  3 of  OH  was 

displaced  by  phosphate  while  OH  of  the  type  1 and  2 are 
unreactive.  The  types  1 and  3 of  OH,  they  considered  as  the 
structural  OH  v;ith  three  vibrations  at  3200,  890,  and 
790  cm 
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Fig.  14.  The  001  face  of  goethite  lattice.  (After  Bragg 
and  Claringbul  1965). 


The  formation  of  binuclear  bridging  resulted  from  the 
displacement  of  two  adjacent  type  3 hydroxyls 
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Fig.  15.  The  001  face  of  phosphated  goethite. 

Phosphate  Adsorption  and  Fesorption  Studies 
Effects  of  the  Initial  Concentration 


The  phosphate  adsorption  isotherm  is  obtained  by  plot- 
ting the  amount  of  phosphorus  adsorbed  (yg  P/g  of  adsorbent) 
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against  the  equilibrium  concentration  of  phosphorus  in 
solution.  The  shape  of  each  isotherm  curve  is  altered  by 
the  relative  amount  of  phosphorus  adsorbed  at  each  equilib- 
rium concentration  and  reaction  time . For  a particular  time 
of  adsorption  reaction,  the  amount  of  phosphorus  adsorbed 
increases  with  the  initial  concentration  but  the  percentage 
of  phosphate  sorbed  decreases. 

For  the  goethite-P  solution  system,  the  shapes  of  the 
phosphate  adsorption  isotherm  curves  were  affected  by  the 
relative  amount  of  phosphorus  adsorbed  at  each  equilibrium 
concentration  (Fig.  16).  Each  curve  was  composed  of  three 
parts:  one  at  low  solute  concentration  (<  0.25  yg  P/ml),  a 

second  part  at  higher  solute  concentration  (0.25  to  1.5  yg 
P/ml)  where  the  isotherm  becomes  convex,  and  a third  linear 
part  at  higher  solute  concentration  (>  1.5  yg  P/m.l). 

However,  in  the  soil-solution  system,  phosphate  adsorp- 
tive capacity  depended  on  soil  characteristics  which  affected 
the  shape  of  the  isotherm  curves.  The  Kenya  soil  and  Georgia 
soil  phosphorus  sorption  curves  apparently  had  two  portions, 
the  first  part  which  is  a curved  portion  (<  5 yg  P/ml)  and  a 
second  or  linear  portion  (>  5 yg  P/ml).  In  Fig.  17,  the 
apparent  lack  of  an  initial  linear  portion  for  the  curve 
found  for  Kenya  soil  suggests  need  for  more  adsorption  data  at 
low  equilibrium  concentration  (<  3 yg  P/ml).  For  the  Georgia 
soil,  (Fig.  18),  the  adsorption  maximum  is  low  compared  to 
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Fig.  17.  Adsorption  isotherm  for  the  Kenya  soil  after  24 
hours  of  reaction  time. 
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Fig.  18.  Adsorption  isotherm  for  the  Georgia  soil  after  a reaction  time  of  24  hours. 
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that  for  the  Kenya  soil,  which  means  that  for  the  same 
initial  concentration  in  solution  more  phosphate  ions  are 
present  in  the  Georgia  soil  than  in  the  Kenya  soil.  For  the 
Colorado  soil,  (Fig.  19),  there  is  a slight  change  in  slope 
for  the  initial  and  final  linear  party  of  the  absorption 
isotherm.  For  the  initial  range  of  equilibrium  concentration 
ranging  from  0 to  25  yg  P/ml,  the  amount  adsorbed  increased 
linearly  as  the  equilibrium  concentration  increases.  This 
relationship  was  20  yg  P adsorped  per  gram  of  soil  for  each 
yg  P/ml.  Muljadi  et  al . (1966)  noted  that  their  isotherm 

curves  were  also  linear  initially  with  a curved  transition 
to  the  second  linear  portion.  They  ascribed  the  initial  and 
curved  portions  to  phosphate  exchange  with  OH  of  Al(OH) 
located  on  the  clay  edge  surfaces.  They  do  not  give  a clear 
explanation  of  the  mechanism  of  adsorption  reaction  for  the 
third  part  of  the  curve  but  postulated  that  the  final  linear- 
ity of  the  isotherm  indicated  that  the  number  of  adsorption 
sites  remained  constant  even  though  the  amount  of  phosphate 
adsorbed  increased. 

The  time  required  for  the  reaction  to  reach  equilibrium 
decreased  as  the  initial  concentration  decreased,  (Fig.  20). 
If  the  initial  concentration  was  less  than  or  equal  to 
5 yg  P/ml,  the  steady  state  of  reaction  was  reached  in  less 
than  10  hours.  Tv;enty  hours  of  reaction  were  required  to 
approach  the  equilibrium  state  with  an  initial  concentration 
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Fig.  20.  Effects  of  the  initial  P concentration  on  the  phosphate  adsorption  by 
goethite-solution  (1  g/1000  ml). 
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greater  than  or  equal  to  5 yg  P/ml.  The  above  observations 
showed  that  low  amounts  of  phosphorus  were  almost  instanta- 
neously adsorbed  onto  the  synthetic  goethite.  Because  the 
initial  phosphate  potential  on  the  solid  phase  was  low  or 
nil  relative  to  the  phosphate  potential  in  solution,  the 
phosphate  flux  from  the  solution  was  therefore  high.  This 
movement  of  phosphate  to  the  solid  surface  continued  for  a 
longer  period  of  time  if  the  initial  concentration  (or 
phosphate  potential)  in  solution  was  high. 

Effects  of  the  Supporting  Electrolyte 
Type  of  Electrolytes 

The  adsorption  isotherm  of  phosphate  on  geothite  was 
examined  using  molar  concentration  of  the  electrolytes, 

NaCl,  CaCl^ , and  AlCl^ , as  given  in  Fig.  21.  Where  the 
initial  concentration  was  greater  than  0.5  yg  P/ml,  the 
salts  (electrolyte)  gave  a significant  effect  on  the  amount 
of  phosphorus  adsorbed.  The  adsorption  at  each  concentration 
decreased  in  the  order  IM  A1C12>1M  CaCl2>lM  NaCl.  It  was 
noted  that  the  adsorption  of  phosphate  using  water  as  support- 
ing medium  gave  the  same  adsorption  isotherm  as  that  using 
1 M NaCl  solution.  For  equilibrium  concentration  greater 

than  1 ug  P/ml,  the  amount  of  P sorbed  in  1 K AlCl  and 

0 

1 M CaCl^  was  greater  than  that  sorbed  in  1 H NaCl  by  a 
factor  of  1.6  and  1.5,  respectively. 


H AICI3 
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Equilibrium  Concentration  (qg  P/rrd) 

Fig.  21.  Effects  of  the  type  of  supporting  electrolyte  on  P adsorption  on 
goethite-suspension  (1  g/1000  ml). 


79 


Table  2.  Type  of  salt  effects  on  the  P adsorption  on 
goethite  at  5 pg  P/ml. 


Electrolyte 

P adsorption  (pg  P/g) 

1 M NaCl 

3780 

1 M -CaCl^ 

5600 

1 M A1C1„ 

5900 

— 3 

Ryden  and  Syers  (1975)  reported  that,  where  some  soils 

have  a final  concentration  above  0.1  pg  P/ml,  the  P sorption 
_ 2 

in  10  M Ca  solution  was  1.5  to  2.5  greater  than  the  sorp- 
tion from  water.  From  the  isotherms  for  P sorption  by 
goethite,  the  data  were  arranged  as  the  plots  of  a/b  versus 
a,  where  a is  the  equilibrium  concentration,  b is  the  amount 
of  phosphorus  adsorbed  per  unit  weight  of  adsorbent.  The 
linearity  of  the  plots  (Fig.  24)  confirmed  the  Langmuir  type 
of  adsorbent.  When  using  Eq . (14c)  the  adsorption  maximum 
(s)  and  the  adsorption  energy  constant  (k)  are  calculated  as 


shown  in  Table  3. 
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Equilibrium  Concentration  (pg  P/ml) 


Fig.  24.  The  transformed  Langmuir  equations  for  phosphate 

sorption  by  goethite  as  affected  by  reaction  times, 
a is  equilibrium  concentration  and  b amount  adsorbed. 
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Table  3.  Effects  of  three  electrolyte  salts  on  the  phos- 
phate sorption  maximum  and  sorption  energy  con- 
stant for  goethite. 


Electrolyte 

Adsorption  Maximum 

Sorption  Energy  Constant 

(yg  P/g) 

(ml/yg  P) 

1 M NaCl 

4300 

3 . 8 

1 M CaCl^ 

6900 

2.3 

1 M AlCl^ 

12700 

1.1 

The  increase  of  phosphate  sorption  maximum,  (Table  3), 
was  accompanied  by  a decrease  in  the  apparent  sorption 
energy  constant.  The  increase  in  phosphate  sorption  may 
be  due  the  increase  in  the  cation  charge  of  the  electrolyte 
favored  adsorption.  A probable  way  in  which  the  cation 
charge  enters  into  the  phosphate  adsorption  reaction  was 
through  cation  bridging: 


Fe  0 


so  that  where  M is  Na,  then  n = 0 , or  if  M is  Ca  then  n = 1, 
and  if  M is  A1  then  n = 2.  As  the  charge  on  the  cation 
increased,  there  v:as  a greater  attraction  between  the  cation 
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on  the  goethite  surface  and  the  phosphate  ion.  The  work 
required  to  bring  the  phosphate  ion  to  the  cation  evidently 
decreased  as  the  valence  increased.  The  increase  in  metal 
valence  favored  a higher  probability  that  the  phosphate  was 
maintained  in  the  goethite-P  solution  interface,  so  that  the 
energy  for  adsorption  decreases. 

Both  aluminum  and  calcium  ions  also  probably  reacted 
with  phosphate  in  solution  to  form  new  phases 
(precipitation),  inducing  thereby  a decrease  of  phosphate  in 
solution.  Som.e  of  the  compounds  which  might  precipitate  in 
solution  depended  both  on  phosphate  concentration  (molarity) 
and  pH.  These  systems  can  be  written  as  follows: 

For  CaHPO^ , the  pH  and  H2P0^  relation  is 

pH2P0^^  = pH  - 3.14  (86) 

For  A1 ( OH ) ^H^PO^ , the  pH  and  H^PO^  relation  is 

pH^FO^^  = pH  + 10.7  (87) 

For  Fe  (OH ) ^H^PO^ , the  pH  and  H^FO^^  relation  is 

pH^PO^^  = pH  - 10.9  (88) 

The  plots  of  pH^PO^^  versus  pH  gave  the  solubility  diagrams 
of  the  compounds  as  illustrated  by  Lindsay  and  Moreno 
(1960).  At  any  point  (pH,  pH^FO^)  above  the  line  for  a 
selected  compound,  precipitation  is  expected  while  below  the 
lines  dissolution  of  the  corresponding  solid  phase  will 
occur.  In  our  study,  the  pH  ranged  from  6 to  7 so  that  with 
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a phosphate  concentration  of  5 yg  P/ml  or  1.72x10“^  moles  of 
H^PO^  per  1 liter,  we  would  have  pH2P0^^  = -log  H^PO^^  = 3.8; 
hence,  no  important  amount  of  precipitation  of  the  above 
solid  phase  was  expected.  Other  complex  compounds  involving 
combination  of  Ca  and  Fe  phosphates  or  A1  and  Fe  phosphates 
could  precipitate  near  or  on  colloid  surfaces  of  the  goethite 
or  the  soils. 

Effects  of  Electrolyte  Concentration 

Salt  concentration  had  a marked  effect  on  the  amount  of 

phosphate  sorbed  during  reactions  at  various  periods.  The 

time  at  which  equilibrium  was  reached  was  evidently  not 

affected  by  increasing  the  salt  concentration  (Fig.  22). 

When  the  initial  concentration  is  5 yg  P/ml,  at  equilibrium 

state  the  phosphate  sorption  in  1 N CaCl2  and  0.01  N CaCl2 

is  1.3  and  1.2,  respectively,  greater  than  that  found  for 

water  system  without  salt.  Van  Olphen  (1977)  pointed  out 

that  increasing  the  electrolyte  concentration  not  only 

caused  compression  of  the  diffuse  part  of  the  double  layer 

but  also  some  ions,  as  counter  ions,  shift  from  the  diffuse 

layer  to  the  Stern  layer.  As  a result,  the  ion  concentration 

of  the  diffuse  layer  decreased  and  more  adsorption  took 

place.  This  concept  was  supported  by  Ryden  and  Syers  (1975) 

who  found  that  equilibrium  phosphate  concentration  remained 

- 2 -1 

the  same  in  both  10  M Ca  and  10  M Na  systems.  Donnan 
theory  could  be  used  to  explain  the  above  phenomenon.  The 
inner  solution  possibly  ranged  from  the  solid  surface  to  the 
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Fig.  22.  Effect  of  the  supporting  electrolyte  concentration 
on  the  kinetics  of  F adsorption  by  the  goethite- 
solution  system  (1  g/1000  ml). 
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Table  4.  Phosphate  desorption  from  goethite  by  different 
anions . 


Anions 

(O.IN) 

Cl 

2- 

SO4 

no; 

cio; 

OH 

H^O 

% Adsorbed  P 
desorbed 

0 . 3 

1.4 

1 — 1 
1 — 1 

CO 
1 — t 

cn 

CO 

2.5 

Table  5.  Effects  of  reaction  time  on  phosphate  adsorption 
maximum  and  sorption  energy  constant  for  the 
goethite  system. 


Reaction  time 
( hours ) 

Adsorption  Maximum 
(ug  P/g) 

Sorption  Energy 
Constant 
(ml/yg  P) 

1/10 

1340 

0.6  6 

1/2 

2280 

0.84 

8 

3850 

2.40 

18  to  76 

4460 

4.60 

Table  6.  The  logarithm  of  (a)  equilibrium  P concentration 
as  a function  of  the  logarithm  of  (b)  amount  of 
P adsorbed. 


log  a 

-.45  .15  3.07  5.54 


log  b 


3 . 40 


3.55 


3.60 


3.65 
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imaginary  limit  of  the  physically  adsorbed  ions  and  the 
outer  solution  containing  all  free  ions.  This  could  be 
expressed  as: 

= (M)J/''(H2P0^^).  = k (89) 

where  (M)^  and  (M)^  are  the  concentrations  of  the  cation  M 
of  valence  n in  the  outer  solution  and  inner  solution  respec- 
tively, terms  (H^PO^^)^  and  (H2P0^)^  are  the  phosphate  concen- 
trations in  the  outer  and  inner  solutions,  respectively,  and 
k is  a constant.  From  the  above  equation,  it  was  evident 
that  as  salt  concentration  (M)^  increased,  the  (H2pO^^)^ 
concentration  must  decrease.  As  a result,  phosphate  may  be 
adsorbed  through  cation  bridging  or  precipitation  as  cation- 
phosphates. 

Phosphate  Desorption 

Effects  of  phosphate  free  solutions  of  1 M AlCl^,  and 
CaCl2  and  NaCl  are  shown  in  Fig.  23.  There  was  a significant 
effect  of  the  type  of  salt  on  the  extent  of  phosphate  desorp- 
tion. The  magnitude  of  differences  due  to  salt  sources 
increased  as  the  equilibrium  phosphate  concentration  in- 
creased. For  a low  equilibrium  concentration  (1  yg  P/ml), 
the  desorption  with  1 M AlCl^  solution  released  550  yg  P/g 
of  goethite.  This  was  3.1  and  1.8  times  greater  than  similar 
desorption  by  1 K CaCl2  and  1 K MaCl,  respectively.  The 
higher  desorption  of  phosphate  by  1 M AlCl^  solution  agreed 
with  the  apparent  lower  potential  binding  energy  constant  of 


Desorption  ( pg  p/g  Goethite) 
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rquilibrium  Concentrntion  (Ug  P/ml) 

Fig.  23.  Phosphate  desorption  in  10  hours  from  goethite 
by  three  salt  solutions  at  various  equilibrium 
F concentration. 
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phosphate  in  goethite-P  solution  in  the  presence  of  this 
salt.  The  increase  in  phosphate  adsorption  due  to  increase 
of  cation  valence  was  also  accompanied  by  an  increase  in  the 
desorbed  fraction  which  was  weakly  bound  (Fig.  23).  This 
would  suggest  that  the  additional  phosphate  adsorption  in 
presence  of  multivalent  cations  may  not  only  be  due  to 
direct  surface  binding  but  also  to  other  bonding  by  some 
form  of  cation  bridging  and/or  multilayer  phosphate  adsorp- 
tion. These  latter  forms  were  easily  replaced  during  P 
desorption  with  further  fresh  salt  solutions. 

Time  of  Reaction  Effects 
Adsorption  Pxeaction 

From  the  decrease  in  the  initial  phosphate  concentration, 
the  phosphate  sorption  increased  with  time  and  remained 
constant  when  the  reaction  time  was  greater  than  18  hours 
for  goethite-P  solution  system  (Fig.  25).  The  adsorption 
made  according  to  the  transformed  Langmuir  equation,  Eq. 

(14c),  was  linear.  This  suggest  that  the  surface  of  the 
synthetic  goethite  was  composed  of  an  homogeneous  population 
of  sites  for  the  adsorption  of  phosphate.  In  Fig.  24,  it 
can  be  seen  that  both  the  slope  and  intercept  of  the  sorption 
isotherms  decreased  as  the  reaction  time  increased  up  to  18 
hours.  In  Table  5 it  v;as  shown  that  both  the  calculated 
adsorption  miaximum  and  the  potential  binding  energy  constant 
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Fig.  25.  Change  in  the  equilibrium  P concentration  and  P sorption  by  goethite  as 
affected  by  the  reaction  time. 
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Fig.  26.  Change  in  the  equilibrium  P concentration  for  three  soils  as  affected  by 
the  reaction  times. 
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increased  as  the  reaction  time  increased.  Changes  in  this 
potential  binding  energy  constant  increased  as  the  reaction 
time  increased.  Changes  in  this  potential  binding  energy 
might  be  explained  by  the  fact  that  as  the  surface  becomes 
covered  by  phosphate,  further  adsorption  of  other  phosphate 
ions  occurred  on  the  remaining  free  surface  which  required 
mueh  more  energy  for  bonding  than  on  the  first  monolayer 
(Clark  1970).  It  should  be  kept  in  mind  that  the  Langmuir 
equation  assumed  that  there  was  only  a monolayer  of  adsorbed 
P and  that  there  was  no  lateral  interaction  between  the 
adsorbed  ions.  As  a result,  the  application  of  the  Langmuir 
equation  to  experimental  data  would  require  that  the  reaction 
time  be  specified  since  the  nature  of  the  adsorbed  layer 
would  not  be  defined  with  time. 

The  increase  of  energy  constant  as  the  adsorption 
maximum  increases  was  not  in  disagreement  with  studies  by 
Aharoni  and  Ungarish  (1977)  who  stated  that  adsorption  took 
place  in  the  initial  stages  on  the  sites  with  lowest  energy 
at  a particular  time.  When  comparing  materials  of  different 
mineralogical  composition,  Ryden  and  Syers  (1975)  found  that 
consistent  trend  for  increase  in  the  energy  constant  with 
the  adsorption  maximum  confirmed  the  assumption  that  phos- 
phate was  usually  adsorbed  onto  all  available  low  energy 
sites  before  the  high  energy  sites  were  filled.  In  the 
present  studies,  using  the  P sorption  data  for  goethite-P 
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solution  when  the  reaction  times  exceeded  18  hours,  further 
treatment  of  the  data  was  made  using  the  transformed  Freund- 
lich  relationship  given  in  Eq . (10).  From  Table  7,  the  plot 
of  log  a versus  log  b was  close  to  linearity  as  the  following 
regression  equation: 

log  b = 0.209  log  a + 3.50  (90) 

Therefore,  the  free  energy  of  adsorption  (expressed  as  the 
energy  difference  between  the  adsorbed  and  free  state  of 
phosphate)  decreased  exponentially  as  the  binding  energy  and 
surface  coverage  increased. 

From  the  proposed  kinetic  model  given  in  Eq.  (57)  and  Eq . 
(59b),  an  attempt  was  made  to  predict  the  change  with  respect 
to  time  that  occurred  for  forms  of  phosphate  logically  in  the 
system.  These  were  (A)  free  ions  such  as  H^PO^,  (B)  the 
physically  and  reversible  chemically  adsorbed  phosphate,  and 
(C)  the  irreversible  chemically  adsorbed  phosphate.  A 
transformation  of  Eq . (57a)  can  be  written  as 

In  a = w^(t  - t^)  + k (91) 

where  a is  the  equilibrium  concentration,  w^  and  k are  con- 
stant, t is  the  reaction  time,  and  t^  is  the  equilibrium 
time  equal  to  18  hours.  According  to  the  sorption,  data  were 
plotted  using  Eq . (91)  and  are  given  in  Fig.  30.  It  can  be 

seen  that  the  phosphate  adsorption  reaction  on  gpethite 
followed  two  steps  before  the  equilibrium  was  reached.  The 
first  step  was  finished  after  3 hours  of  reaction  tim.e  (t^). 


In  (Equilibrium  Concentration  [pg  P/ml]) 
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Fig.  30.  Logarithmic  plot  of  equilibrium  P concentration 

change  with  time  (t)  relative  to  equilibrium  at  18 
hours  (t^)  for  phosphated  goethite.  Data  were 
obtained  using  Eq . (91). 
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and  the  second  step  after  3 to  18  hours  (t^).  From  these 
data  it  was  apparent  that  Eq.  (57a)  should  be  split  into  two 
components  with  respect  to  time  after  initiating  the  experi- 
ment so  that 


a,  = a e 
1 o 


(92a) 


and  for  the  range  t^  to  t^ 


a , e 
ol 


w,,  (t  - t ) 
lb  o 


(92b) 


so  that  a^  is  the  initial  P concentration  in  solution  for 
t ^ and  a2  is  the  solution  P concentration  at  longer 
reaction  time  between  t^^  and  t^.  In  the  goethite-P  solution 
system,  there  was  an  initial  fast  reaction  so  that 
w la  = hour  ^ and  a second  step  of  the  reaction  which 

was  slower  with  = 0.025  hour  This  would  also  apply 

for  the  fast  and  slow  adsorption  reactions  observed  by 
Barrow  and  Shaw  (1975). 


Desorption  of  Phosphate 

After  the  phosphate  adsorption  was  conducted  for  dif- 
ferent periods  of  time,  the  desorption  curve  shown  in  Fig. 
27  for  the  Kenya  soil  indicated  that  the  phosphate  desorp- 
tion in  0.01  N CaCl^  was  almost  linear  with  respect  to  P 
concentration  in  the  equilibrium  solution  and  that  the 
adsorption  reaction  tim.e  affected  the  desorption  process. 
With  an  initial  concentration  of  20  ug  P/ml  after  23  hours 


Fig.  27.  Effect  of  reaction  time  and  equilibrium  F concentra- 
tion on  P desorbed  by  0.01  b CaCl^  from  Kenya  soil. 
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of  adsorption  reactions  in  0.01  N CaCl2  solution  desorption 
was  80  ug  P/g  for  Kenya  soil.  Similarly  at  each  equilibrium 
concentration,  the  desorption  of  phosphate  increased  for  the 
Georgia  soil.  Fig.  28  also  showed  more  desorption  of  P from 
the  Colorado  soil  after  23  hours  of  P adsorption  than  that 
after  27  hours.  For  the  Colorado  soil,  the  phosphate 
desorption  data  indicate  that  if  the  equilibrium  concentra- 
tion was  greater  than  5 pg  P/ml,  there  was  a divergence  in 
phosphate  release  as  the  adsorption  reaction  time  increased. 
In  Fig.  29  extraction  with  0.5  N NH^^F  at  pH3 . 8 resulted  in  a 
significant  increase  in  phosphate  release,  presumably  from 
that  in  aluminum  form.  Apparently,  more  phosphate  was  held 
as  an  aluminum  phosphate  after  72  hours  of  adsorption  than 
after  21  hours,  a negligible  amount  being  recovered  after  4 
hours  of  adsorption  time.  This  might  explain  the  lower 
desorption  of  phosphate  by  0.01  N CaCl2  at  longer  adsorption 
time . 

In  the  goethite-P  solution  system,  (Fig.  25),  phosphate 
adsorption  at  first  increased  and  then  remained  constant  as 
the  reaction  time  was  increased.  This  indicated  that  there 
is  first  a predominant  physically  adsorbed  P followed  by  a 
conversion  to  chemically  adsorbed  P;  this  is  supported  by 
the  increase  in  P sorption  energy  constant  (Table  5).  The 
free  phosphate  will  replace  some  of  the  physically  adsorbed 
phosphate  until  a stable  layer  will  be  established  and  the 
system  is  then  at  equilibrium. 


Desorption  (pg  p/g  Soil) 
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Fig.  28,  Effect  of  reaction  time  and  equilibriumi  F concentra- 
tion on  P desorped  by  0.01  N CaCl„  from  the  Georgia 
and  Colorado  soils. 
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Fig.  29.  Effect  on  P desorbed  by  0.5  N NH^F  as  affected  by 
the  adsorption  reaction  time  and  equilibrium  P 
concentration . 
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Effects  of  pH  on  P Adsorption  and  Desorption 
Adsorption 

The  adsorption  of  phosphates  both  on  goethite  and  soils 
induced  an  increase  in  the  equilibrium  solution  pH,  (Table 
7).  The  calculated  pH  change  was  0.75xl0”^pH  units/pg  P/g 
of  goethite  in  the  soil-P  solution  systems;  it  was  3.75X10" 
2.14X10  and  0.75X10  ^pH  units  per  pg  P/g  soil  for  the 
Georgia,  Kenya,  and  Colorado  soils,  respectively.  In  the 
goethite-P  solution  studies,  phosphate  adsorption  increased 
with  decreasing  pH  of  the  equilibrium  solution  (Fig.  31). 
When  the  initial  P solution  was  5 pg  P/ml,  the  amount  of 
phosphate  adsorbed  decreased  almost  linearly  with  pH  of  the 
equilibrium  solution  was  increased  according  to  the  equation 
(pg  P/g)  = -446  pH  + 5680  (93) 

where  P is  the  amount  of  P absorbed  on  goethite. 

The  effects  of  pH  on  the  prevalence  of  the  orthophos- 
phate species  in  solution  have  been  investigated.  From 
Lindsay  and  Vlek  (1977),  each  of  the  following  phosphate 
species  constituted  more  than  50%  of  the  total  P at  the 
following  respective  pH  range:  i)  H^PO^  for  pH  2 or  less, 

ii)  H^PO^^  for  pH  2 to  7 , iii)  HPO^  for  pH  7 to  12,  and 
3 - 

iv)  P0[^  for  pH  greater  than  12.  Such  changes  in  phosphate 
ions  with  pH  further  complicated  the  interpretation  of 
phosphate  adsorption.  As  the  pH  increased,  the  goethite  had 
an  increasing  net  negative  charge  (Fq.  64).  Since  the 
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Fig.  31.  Effect  of  change  in  equilibrium  solution  pH  on  P 
adsorption  and  0.01  N CaCl_  desorption  of  P in 
goethite  system.  Initial  P concentration  is 
5 yg/m.l. 
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Table  7. 

Effects  of 
solution  pH 

P ' adsorption  on 
after  24  hours 

the  equilibrium 
of  adsorption  reaction 

Adsorbent 

P added 
(pg  P/g) 

P adsorbed 
(yg  P/g) 

pH 

Initial 

Final 

Goethite 

5 

3500 

7.10 

7 .30 

10 

4300 

6 .60 

7 . 00 

Georgia 

5 

179 

5.10 

5 .55 

10 

202 

5.10 

5.47 

Colorado 

5 

198 

6 .80 

7 .03 

10 

374 

6 .57 

6 . 69 

Kenya 

5 

304 

5.50 

6 . 06 

10 

490 

5 .40 

5.95 
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predominant  orthophosphate  ions  are  in  the  increasing  order 
2-  3- 

HPO^  , PO^^  with  increasing  pH,  from  the  electrostatic 
point  of  view,  there  was  an  increasing  repulsion  between  the 
goethite  surface  and  the  orthophosphate  ions.  This  resulted 
in  less  phosphate  adsorbed  with  an  increase  in  pH.  Hingston 
et  al.  (1967)  found  that  the  maximum  adsorption  of  goethite 
was  approached  when  the  pH  was  less  than  4.  At  such  low  pH, 
goethite  would  have  a net  positive  surface  and  phosphate  ion 
would  be  negatively  charged.  As  a result,  minimum  repulsion 
or  maximum  attraction  between  goethite  surface  and  phosphate 
ions  should  occur  at  this  acidity. 

Desorption 

The  reversibility  of  the  adsorbed  phosphates  from 
solution  at  different  pH  condition  is  shown  in  Fig.  31.  The 
pH  at  which  the  phosphate  adsorption  took  place  had  an 
important  effect  on  the  reversibility  of  phosphate  adsorption. 
As  the  pH  of  the  equilibrium  solution  (adsorption  process) 
increased,  there  was  first  a decrease  in  the  quantity  of 
phosphate  desorbed  until  a minimum  was  reached  in  the  range 
pH  5 to  pH  6 followed  by  descreased  desorption  as  pH  was 
further  increased.  At  less  than  pH  5,  both  phosphate  adsorp- 
tion and  desorption  were  greately  increased.  This  indicated 
that,  as  the  pH  decreased,  there  were  some  additional  posi- 
tive sites  formed  on  the  goethite  surface.  One  probability 
was  the  shift  of  hydronium  ion  to  the  hydroxyl  group  of 
Fe  - OH  as  shov;n  below 
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OH 

H - ' 

Fe  - OjJ+  . . . 0 - P = 0 (94) 

OH 

resulting  in  more  phosphate  adsorption.  However,  when  the 
pH  was  greater  than  6,  there  was  an  increase  in  the  hydroxyl 
ion  concentration  which  has  strong  affinity  for  any  positive 
site.  Consequently,  the  existence  of  competition  between 
OH  and  the  orthophosphate  ions  for  the  same  type  of  sites 
resulted  in  lower  phosphate  adsorption.  Hingston  et  al. 
(1968)  noticed  that  anions  were  desorbed  by  competitive 
effect  when  the  competitor  can  occupy  sites  in  addition  to 
those  already  occupied  by  other  anions,  hence  the  increase 
in  the  negative  charge  of  the  surface.  This  increase  in 
negative  charge  produced  by  the  competitor  allowed  phosphate 
desorption  by  surface  hydroxyls.  It  was  known  that  the 
increase  in  OH  concentration  favored  the  formation  of 
surface  negative  sites,  which  contribute  to  the  weak  phos- 
phate binding  resulting  in  greater  P desorption. 

The  competitive  capacity  of  hydroxyl  for  phosphate 
sites  was  further  confirmed  from  data  shown  in  Fig.  32.  The 
pH  effect  on  the  equilibrium  P concentration  after  various 
adsorption  times  was  m.arkedly  different  for  pH2  and  pHlO. 

At  pH2,  the  desorption  of  phosphate  increased  initially  and 
remained  constant  after  5 hours  of  reation  tim.e.  However  at 
pHlO  the  phosphate  desorption  increased  sharply  initially 


Desorption  (ug  P/g  Goethite) 
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Fig.  32.  Effect  of  phosphate  adsorption  time  on  phosphate 
desorption  by  0.01  N CaCl2  at  pH  2 and  pH  10. 
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and  thereafter  decreased  slowly  and  linearily  as  the  reaction 
time  increases.  This  shows  that  at  pfIlO  the  hydroxyls  ions 
continued  the  phosphate  desorption  as  the  adsorption  reaction 
time  was  increased.  This  apparent  increase  in  efficiency 
for  OH  replacement  of  phosphate  might  also  have  included 
calcium  phosphate  precipitation  as  well  as  the  fact  that 
much  of  the  phosphate  ions  were  probably  in  trivalent  state. 
Obviously,  ageing  of  the  phosphated  goethite  negated  phos- 
phate desorption  by  OH  replacement. 

Effects  of  pH  and  Time 

In  order  to  define  the  response  of  phosphates  desorption 
with  time  and  pH,  the  desorption  was  conducted  at  both 
different  periods  of  time  and  pH.  From  the  data  shown  in 
Fig.  33,  at  both  6 and  12  hours  of  desorption  time,  phosphate 
desorption  followed  the  same  pattern  over  the  same  pH  range. 
However,  after  20  minutes,  only  at  pH  10  was  substantial 
amounts  of  phosphate  desorbed  (600  yg  P/g  goethite)  which  is 
about  1/3  of  the  amount  that  was  obtained  after  12  hours  of 
desorption  at  this  pH.  Generally,  it  can  be  said  that  the 
efficiency  of  phosphate  desorption  increased  as  the  pH  of 
the  desorbing  solution  was  increased. 

The  effects  of  the  electrolyte  concentration  and  pH 
interaction  on  phosphate  desorption  are  shown  in  Fig.  34. 

The  phosphate  desorption  at  any  pH  increased  when  supporting 
electrolyte  was  present  compared  to  use  of  water  only  for 
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Solution  pH 


Fig.  33.  Effects  of  solution  pH  on  the  amount  of  F desorped 
from  goethite. 
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Fig.  34.  Effects  of  pH  and  concentration  of  supporting  elec- 
trolyte on  phosphate  desorption. 
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phosphate  desorption.  In  the  range  pH  2 to  pH  12,  either 
1 N CaCl2  or  0.01  N CaCl2  gave  similar  amounts  of  desorption 
phosphate  at  each  pH,  and  each  is  greater  than  that  obtained 
from  water.  At  any  pH,  the  phosphate  desorption  in  0.01  N 
CaCl2  was  approximately  200  yg  P/g  goethite  greater  than  the 
amount  desorbed  in  water. 

Effects  of  P Adsorption  on  Surface  Charge  of  Goethite 
In  goethite-solution  and  soil  solution  system,  it  was 
assumed  that  H and  OH  ions  are  potential  determining  ions. 
Potentiometric  titration  of  goethite  was  made  to  determine 
the  zero  point  of  charge  (ZPC)  from  the  distribution  of  the 
net  electric  charges  with  varying  pH  and  electrolyte  concen- 
trations (Fig.  35).  The  ZPC  was  at  the  inflexion  point  of 
the  titration  curves  at  different  electrolyte  concentrations. 
Parfitt  and  Atkinson  (1976)  reported  a a decrease  of  the  pH 
for  ZPC  from  8.1  to  5.1  as  a result  of  the  adsorption  of  100 
ymole  of  NaH2P0^  on  goethite.  According  to  Mekaru  and 
Uehara  (1972),  each  millimole  of  P sorbed  increased  the  CEC 
by  approximately  0.8  meq  at  pH  7 on  some  ferruginous  tropical 
soils.  In  this  study,  when  using  synthetic  goethite,  ZPC 
was  at  pH  5.8  for  goethite  (shown  in  Fig.  36),  compared  to 
ZPC  at  pH  5.2  for  phosphated  goethite  (as  a result  of  3000 
yg  P adsorbed  per  g goethite).  From  Eq . (66)  it  was  deduced 
from  Gouy-Chapman  theory  that  there  should  be  an  increase  in 
negativity  as  a result  of  phosphate  adsorption.  In  the 


Adsorbed  (meq  per  100  g Goethite ) 
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Fig.  35.  Fctentiometric  titration  of  goethite.  Note 
point  of  charge  occurs  at  pH  5.8. 


zero 


Sdsorbed  (meq  per  100  g Goethite) 
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Fig.  36.  Fotentiometric  titration  of  phosphated  goethite. 
Not  zero  point  of  charge  occurs  at  pFI  5.2. 
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present  study,  the  goethite  was  freed  of  electrolyte  by 
dialysis  prior  to  ZPC  determinations.  This  process  may  have 
had  some  effect  on  the  nature  of  the  goethite  surface  since 
loosely  bound  ions  were  removed  in  the  process.  In  a soil 
system,  all  surfaces  have  various  cations  and  anions  present 
on  the  clay  surfaces  which  affect  both  phosphate  adsorption- 
desorption  and  ion  activities  resulting  after  phosphate 
adsorption.  From  the  present  studies,  both  pH  and  electro- 
lyte changes  greatly  affected  the  degree  of  phosphate 
desorption . 


CHAPTER  VI 

SUMMARY  AND  CONCLUSION 


The  low  availability  of  phosphate  ions  in  tropical 
soils  is  of  major  concern  to  all  workers  recognizing  the 
need  to  increase  agricultural  production.  In  this  study, 
conditions  favoring  the  formation  of  iron  oxyhydroxides  such 
as  goethite  were  investigated.  An  attempt  was  made  to 
understand  the  mechanism  through  which  phosphate  may  be 
bound  to  goethite  using  infrared  spectroscopy.  Using 
goethite  and/or  soils  as  supporting  medium,  the  influence  of 
some  factors  on  the  kinetics  of  phosphate  adsorption  and 
desorption  were  studied. 

The  conditions  for  goethite  preparation  were  investi- 
gated. Goethite  was  made  by  mixing  FeCl^  and  NaOH  solutions 
and  allowing  the  precipitates  to  crystallize.  The  identifi- 
cation of  the  presence  of  goethite  in  the  suspension  was 
made  by  infrared  spectroscopy  supported  by  X-ray  diffraction. 
The  infrared  spectra  indicated  that  the  bands  characteristic 
of  goethite  formation  are  at  3200  cm"^  (Fe  - OH  stretching), 
890  cm  ^ and  790  cm  ^ (which  are  both  Fe  -OH  bending  vibra- 
tions). The  appearance  of  2.19  8 and  2.70  ^ peaks  from  the 
X-ray  diffraction  also  was  supporting  evidence  that  most  of 
the  components  resulting  from  mixtures  with  OH/Fe  = 6 were 
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goethite.  Both  the  speed  and  the  degree  of  goethite  crystal- 
lization increased  as  the  OH/Fe  ratio  was  increased.  This 
was  seen  through  the  increase  in  the  above  band  intensities. 
From  this  study,  goethite  was  not  a major  component  in  a 
suspension  of  OH/Fe  = 1.5  which  was  aged  only  for  50  hours, 
as  reported  by  some  workers.  With  OH/Fe  = 1.5  in  the  suspen- 
sion, there  was  need  both  to  age  for  a longer  period  of  time 
and  to  increase  the  temperature  of  the  suspension  to  increase 
the  formation  of  Fe  - OH  bonding. 

The  presence  of  phosphate  at  the  beginning  of  the 
ageing  process  weakened  the  structure  of  goethite.  From  the 
data  found  experimentally,  evidence  of  competition  existed 
between  the  phosphate  and  the  hydroxyl  ions  for  bonding  with 
iron.  The  increase  of  phosphate  concentration  in  the  suspen- 
sion (P/Fe  - 0.032  to  3.2)  induced  more  and  more  the  forma- 
tion of  the  P - 0(Fe)  bonds  as  seen  by  the  1000  cm~^  vibra- 
tions. The  surface  adsorption  of  phosphate  on  goethite  was 
identified  through  the  appearance  of  vibrational  bands'  in 
the  1200  - 1000  cm  ^ region.  The  appearance  of  hydroxyl 
deformation  in  the  same  region  was  also  identified  by  treat- 
of  the  phosphated  goethite  with  D2O.  However,  it  was 
believed  that  the  bands  at  1160,  1085,  and  1030  cm"^  for 
phosphated  goethites  were  characteristic  of  phosphate  binu- 
clear  bridging,  (Fe0)2P00H.  The  existence  of  such  type  of 
binding  is  further  confirmed  by  the  increase  of  negative 
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surface  charge,  or  decrease  of  the  point  of  zero  charge 
after  the  phosphate  adsorption. 

Using  synthetic  goethite,  phosphate  adsorption  and 
desorption  were  affected  by  the  initial  concentration,  the 
supporting  electrolyte,  the  reaction  time,  and  pH  of  the 
system.  The  amount  of  phosphate  adsorbed  was  proportional 
to  the  initial  phosphate  concentration.  In  goethite-solution 
system  (1  g/1000  ml),  maximum  phosphate  adsorption  was 
reached  when  the  initial  solute  concentration  was  at  least 
10  yg  P/ml.  The  calculated  adsorption  maximum,  according  to 
the  Langmuir  equation,  increased  in  each  case  with  time,  the 
electrolyte  concentration,  and  the  valence  of  the  cation. 
Increase  in  phosphate  adsorption  resulted  from  an  increase 
in  salt  concentration.  This  was  also  accompanied  by  a 
substantial  increase  in  phosphate  desorbed  by  water.  The 
increase  in  phosphate  adsorption  due  to  increasing  cation 
valence  was  also  accompanied  by  an  increase  in  phosphate 
desorption  by  the  water  in  the  corresponding  electrolyte 
over  that  for  water  only.  The  addition  of  cations  might 
have  contributed  to  the  phosphate  adsorption  either  through 
cation  bridging  or  to  a multilayer  of  physically  adsorbed 
phosphate.  This  was  supported  by  the  fact  that  an  increase 
of  cation  valence  was  not  accompanied  by  an  increase  of  the 
apparent  sorption  energy  constant.  Only  with  time  was  the 
phosphate  adsorption  reaction  accompanied  by  an  increasing 
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sorption  energy  constant.  Changing  of  the  goethite  suspen- 
sion and  phosphate  systems  from  pH  2 to  pH  12  resulted  in  a 
change  in  phosphate  adsorption  because  there  was  competition 
between  the  hydroxyl  and  phosphate  ions  for  adsorption 
sites.  The  adsorption  decreased  almost  in  a linear  fashion 
as  the  pH  increased.  At  the  same  initial  equilibrium  P 
concentration,  the  phosphate  desorption  was  reduced  when  the 
adsorption  solution  pH  was  less  than  6;  but  when  the  pH  was 
greater  than  6 the  desorption  phosphate  increased.  This 
showed  that  more  phosphate  was  displaced  from  the  chemically 
adsorbed  forms  to  desorbable  forms.  This  phosphate  desorp- 
tion increased  with  time  up  to  6 hours  and  remained  constant 
from  6 to  12  hours. 

Both  Langmuir  and  Freundlich  equations  were  used  to 
describe  the  relationship  between  the  amount  of  phosphorus 
absorbed  and  the  equilibrium  concentration  in  goethite  and 
soils  systems.  Since  those  equations  were  valid  only  at  a 
particular  reaction  time,  a kinetic  model  was  proposed  to 
describe  the  change  of  ions  in  solution  as  the  reaction  time 
increased.  The  model  fractionated  phosphate  absorption 
change  with  time  for  the  physically  adsorbed  and  reversible 
chemically  adsorbed  phosphate  ions  compared  to  that  for  the 
irreversible  chemically  adsorbed  phosphate.  From  the  proposed 
•kinetic  model,  an  expression  for  minimum  energy  for  irrevers- 
ible fixation  was  deduced.  The  experimental  data  showed 
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that  before  the  equilibrium  state  was  reached,  the  adsorp- 
tion reaction  proceeded  through  two  essentially  linear-type 
reactions  with  respect  to  the  reaction  time. 

The  phosphate  adsorption  isotherm  for  soils  indicated 
no  major  differences  from  those  of  other  workers.  The  Kenya 
soil  had  the  highest  adsorption  compared  to  that  of  a Colo- 
rado soil  and  a Georgia  soil.  Phosphate  adsorption  reached 
the  equilibrium  state  after  22  to  24  hours  of  reaction  time. 
For  any  initial  concentration  greater  than  5 yg  P/ml,  desorp- 
tion of  phosphate  by  0.01  N CaCl^  solution  increased  with 
the  absorption  reaction  time  for  goethite  and  for  the  Kenya 
and  Georgia  soils,  but  not  for  the  Colorado  soil.  This 
desorption  was  in  agreement  with  the  fact  that  after  phos- 
phate desorption  by  0.01  N CaCl2,  further  phosphate  extracted 
by  the  0.5  N NH^^F  increased  as  the  adsorption  reaction  time 
increased.  Some  of  the  phosphate  adsorption  at  shorter  time 
periods  was  extractable  by  0.01  N CaCl^,  but  at  longer 
adsorption  reaction  periods,  phosphate  bonding  strength 
increased. 

Some  of  the  factors  determining  the  phenomena  of 
adsorption-desorption  have  been  investigated.  In  order  to 
relate  these  studies  to  phosphate  fertilization  of  tropical 
soils  and  phosphate  adsorption-desorption  reaction,  such 
studies  conducted  in  the  laboratory  should  be  complemented 
by  pot  studies  involving  tropical  plants.  Later,  from  these 
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findings,  field  experiments  would  be  necessary  to  provide 
demonstration  plots  on  advances  in  technology  for  phosphorus 
fertilization  of  tropical  soils. 
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